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Chapter 1 

HSCT Aircraft Emissions Data* 

SUMMARY 

Estimates are given for the emissions from a proposed High Speed Civil Transport (HSCT). 
This advanced technology supersonic aircraft would fly in the lower stratosphere at a speed in the 
range of roughly Mach 1 .6 to Mach 3.2 (470 to 950 m/sec or 920 to 1850 knots). Because it vould 
fly in the stratosphere at an altitude in the range of 15 to 23 km commensurate with its design 
speed, the potential exists for its exhaust effluents to perturb the chemical balance in the upper 
atmosphere. The first step in beginning to determine the nature and magnitude of any chemical 
changes in the atmosphere due to these proposed aircraft is to identify and quantify the chemically 
f important species they emit. 

This chapter summarizes relevant earlier work dating back to the Climatic Impact Assessment 
Program (CIAP) studies of the early 1970’s and current propulsion research efforts at NASA and 
at its High Speed Research Program (HSRP) contractors’ laboratories. Recent work funded by 
HSRP is providing estimates of the chemical composition of an HSCT’s exhaust and these emission 
indices (Els) are presented. Other aircraft emissions that are not due to combustion processes are 
also summarized and these emissions are found fo be much smaller than the exhaust emissions. It is 
anticipated that future advances in propulsion technology, in experimental measurement techniques, 
and in understanding of upper atmospheric chemistry may affect these estimates of the amounts of 
trace exhaust species or their relative importance, and revisions are expected to become necessary 
in the future. 

INTRODUCTION AND BACKGROUND 

The upper atmosphere has been perturbed by anthropogenic chemicals. Most notable of these 
are the chlorofluorocarbons (CFCs), which are implicated in the depletion of stratospheric ozone 
and the creation of the Antarctic ozone hole. Because of the variety of technical advantages 
they offered, these chemically inert species were developed for industrial processes and consumer 
products before it was realized that they were photo-chemically reactive in the upper atmosphere. 
While an understanding of the detailed chemical balance and how it is affected by perturbations 
due to pollutants is still being developed, it is with an awareness of actual anthropogenic effects on 
the stratosphere that the exhaust emissions from the proposed HSCTs are being analyzed. 

The exhaust species under most scrutiny at present are the oxides of nitrogen, NO and NO 2 , 
collectively denoted NO x . The potential for catalytic destruction of ozone by exhaust NO x in the 
stratosphere was recognized 1 in the early 1970’s when the U.S. supersonic transport (SST) was 
being studied. The Climatic Impact Assessment Program (CIAP) was undertaken to consolidate 
and extend existing knowledge of the chemistry, physics, and technology of supersonic flight in the 
stratosphere and the effects of the consequential exhaust emissions. The results and conclusions of 
that program are summarized in its proceedings 2-5 and monographs 6-11 . The state of knowledge 
and technology on the present subject of aircraft emissions as of 1975 is represented in the Second 
Monograph 7 on ‘Propulsion Effluents in the Stratosphere’ . 


* This chapter also is published as Chapter 2 of NASA Reference Publication 1272, “The Atmospheric Effects 
of Stratospheric Aircraft: A First Program Report” 
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Aircraft propulsion technology has advanced in the intervening years such that the reasons the 
U.S. SST was not considered viable can now be addressed with technical improvements 12-14 . The 
evolution of such technically improved engines has actually made NO x reduction more difficult 
as the combustor pressure and temperatures have risen to improve propulsion efficiency. Yet 
reduction of NO x levels below that achieved with SST-era technology appears to be necessary 
to avoid large stratospheric ozone depletion. In order to do this, combustors must be improved, 
requiring the results from current research aimed at controlling local temperatures and equivalence 
ratios throughout the combustion process. A pivotal question concerns the quantitative effects on 
the atmosphere of the reduced levels of exhaust trace species emitted from the advanced technology 
propulsion systems that would power the proposed HSCT. 

The understanding of stratospheric chemistry and transport has also made enormous strides. 
The role of heterogeneous chemistry, in particular, was not apparent when commercial stratospheric 
flight was last considered. An awareness of chemistry on the surface of condensed water and/or 
condensed aqueous solutions has grown out of the need to understand the Antarctic ozone hole. 
Now this heterogeneous chemistry must be accounted for in the global stratospheric chemical 
balance, but also must be assessed for a possible role in the wake of stratospheric aircraft. In the 
aircraft wake, the locally high (relative to ambient concentrations) trace species concentrations 
may encounter condensed water in the condensation (contrails) in the exhaust behind the aircraft. 

Additional data concerning condensation nuclei (CN) emitted by the aircraft are needed to 
understand the possible role of heterogeneous chemistry, both globally and in the wake. These CN 
are essential for the formation of contrails and their number density affects the contrail particle sizes 
and number densities. The chemical nature of the CN surfaces control their ability to condense 
water, newly formed soot typically having a small fraction of its number density as active CN. 

CN in the form of ambient condensed sulfate parti cles, exhaust carbonaceous soot particles, and 
possibly other exhaust particulates are necessary to initiate the condensation process that can occur 
under some conditions in the proposed stratospheric flight paths. Measurements of CN number 
densities and knowledge of their condensation properties are necessary to predict droplet lifetimes 
and settling distances and thus their role in wake chemistry and transport and in global stratospheric 
aerosol loading. 

Once a condensed surface is present, additional chemical species must be accounted for, and 
additional chemical reactions occurring on the condensed phase/vapor interface must be introduced 
in the overall chemical balance. HN0 3 and N 2 0 5 formed by further reactions of NO x , have been 
shown to figure prominently in the heterogeneous chemistry occurring in polar stratospheric clouds 
(PSCs) that drives the Antarctic ozone hole. These species c ould also react with aerosols in the 
contrail if and when concentrations of species and particles are sufficiently high for long enough 
times. 

The list of chemically relevant exhaust species now goes beyond the major combustion products 
of C0 2 and H 2 0 and the trace species NO x . A fuller set of ‘odd- nitrogen’ compounds termed 
NO y , including NO x , NO 3 , N 2 0 5 , HNO 3 , and (although not an exhaust species) C10N0 2 , must 
be considered as aerosol-active species, as well as S0 2 and soot particles. The total unbumed 
hydrocarbons (THC) and CO in the exhaust represent combustion inefficiencies and play a role in 
important stratospheric HO x (OH, H0 2 ) chemistry. 
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The measurement techniques used to measure these species will be discussed in the next section 
followed by a brief summary of how engine emissions are being projected for the next generation of 
low NO x , sustained supersonic propulsion systems. The projected estimates will be discussed in a 
following section and compared with past measurements on related, predecessor engines. Finally, 
a set of estimates for emissions from the aircraft that are not due to the propulsion system \ ill be 
reviewed. 

MEASUREMENT APPROACHES 

Quantitative measurements of exhaust emissions are carried out only in ground-based f cilities. 
In fact, emissions have rarely been measured during aircraft flight, and only qualitative results 
were obtained 15-17 . Ground-based measurements of entire engines must be performed in a test 
facility that can reproduce ambient conditions, including pressure and temperature, if quantities 
appropriate to high altitude flight are to be measured. Such tests are exceedingly expensive 18 
and require considerable preparation and support, so there is strong motivation to perform simpler 
measurements on components of the engine separately in smaller laboratory based facilities. 

Primarily, these simpler measurements consist of reproducing the flow conditions in a trial 
combustor or a simplified version of one and making measurements at its exit. The flow into 
and out of the combustor can be calculated reliably for a chosen engine cycle. The more benign 
conditions encountered during laboratory component testing increase the accuracy of individual 
measurements. 

Most gaseous species are measured using continuous sampling probes that take a small flow of 
the exhaust gases from the component to the measurement instrument, controlling the temperature 
and flow velocity, and thus the delay to measurement. The’ species currently being measured 
routinely include C0 2 , H 2 0, CO, NO, N0 2 , S0 2 ,andthe total unburned hydrocarbons (THC) 19-20 , 
both as an undifferentiated sum of total C as well as individual hydrocarbon species in more 
detailed batch sampling studies 21-23 . Batch samples are analyzed using gas chromatography 
and/or mass spectroscopy to quantify the various hydrocarbon species present, as well as H 2 , 0 2 , 
N 2 , CO, and C0 2 . Particulates in the form of carbonaceous soot have been a concern for decades; 
originally attempts were directed toward reducing the visible smoke leaving with the exhaust. 
These particulate emissions have typically been reported as an SAE ‘smoke number’ (SN) 20 that 
can be related to total particulate mass fairly reliably 24 but, when particle size distributions, number 
densities, and CN-activity are needed,this measurement approach is not adequate. More detailed 
sizing measurements have been performed 25 , but are not routine. 

Chemiluminescence detectors have been used to measure NO as well as NO x . For NO x 
measurements, the N0 2 fraction of the NO x is first converted to NO and the chemiluminescence 
generated when the net NO reacts with ozone is quantified 19 . N0 2 is determined from the difference 
between an NO measurement and the corresponding NO x measurement. NO x measurements have 
conventionally been reported as the total mass of N0 2 plus the mass equivalent of NO ox iized 
to N0 2 . Thus NO x is conventionally reportecTas the N0 2 gram equivalent of the total mea rred 
species 19,26 , independent of the oxidation state (even though it is typically mostly NO near the 
exit of gas turbine engines). This convention, unfortunately, has not always been stated explicitly, 
generating some confusion in interpreting these numbers in the broader community. In this chapter, 
NO, N0 2 , and NO x emissions will all be reported as gram equivalents of N0 2 . Generally, other 
NO y species (e.g. HN0 3 ) have not been measured to determine if significant levels are present. 
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While chemiluminescent N0 X measurements have been and continue to be the standard NO x 
measurement technique, ultraviolet absorption measurements of NO were made as part of the CIAP 
program 27 . A large discrepancy was observed between NO measured using absorption and NO x 
quantified by collecting exhaust with a relatively crude probe and analyzing it with a chemilumi- 
nescent instrument; the UV absorption being up to 4-5 times larger than the chemiluminescent NO x 
measurement 28 . This discrepancy was largest in high temperature and high velocity flows such as 
would be encountered in the supersonic exhaust from any HSCT engine. Improvements in probe 
design and line-of-sight measurement methods may allow these differences to be resolved, but a 
quotable reference study addressing the discrepancies between optical and sampling measurements 
appears to be lacking. Given the need for accurate and reliable emissions estimates required for 
modeling the stratospheric impact of these vehicles, these differences must be understood and an 
optimal measurement strategy adopted. 

Flame ionization has been used to quantify the total unbumed hydrocarbons (THC), while non- 
dispersive infrared techniques are used to measure CO, C0 2 , and water vapor 19 . S0 2 is directly 
related to the sulfur content of the fuel and thus is not controlled by combustor design. When S0 2 
is measured, a non-dispersive ultraviolet technique is employed. The SAE SN is measured for 
particulates emissions by collecting particles as a known volume exhaust sample is drawn through 
a filter and quantifying the change in the filter’s reflectance 20 . 

Radical species, particularly OH, may be chemically active in the exhaust leaving the engine 
nozzle. Although it is not one of the standard exhaust species usually reported, superequilibrium OH 
concentrations were measured in the exhaust of a turbojet for supersonic flight conditions during 
CIAP 27 . The observed superequilibrium ratio was large (> 10) in non-afterburning operation, 
while near unity for afterburning cases. Thus OH concentrations and chemical effects are expected 
to be dependent on engine configuration and cycles (duct-burning versus turbofan/turbojet). If 
an important role in exhaust chemical processing of NO y and CN is established, some reliable 
measurement of OH would be necessary. 

The measurement of a species in the exhaust flow dete rm ines its local concentration or mass 
fraction, but with the variety of combustor designs, engine cycles and bypass ratios being considered 
for these engines, different amounts of air end up in the exhaust for a given fuel flow rate. It has 
been convenient to normalize the emission rate, on a mass basis, of a given species by the mass 
flow rate of fuel. Thus an emission index (El X) is defined as 


ei (Y) x = 


(g/s species X) 
(kg/s fuel) 


( 1 ) 


Emissions expressed as volume fractions require conversion using species molecular masses 19 
in order to arrive at an El. For nearly complete combustion, calculation of an emission index 
effectively normalizes by the rate of enthalpy addition to the air as it passes through the engine. 

The preceding subscript (Y) is not a standard notation but will be used in this chapter to 
indicate that the emission index of species X is being reported as the massequivalent of species Y. 
This notation resolves any possible ambiguities associated with reporting NO x Els and allows the 
reporting convention in terms of N0 2 to be explicitly indicated. THC is also a mixture of compounds 
and its Els are reported in terms of CH 2 , CH a (where a is the fuel atomic hydrogen/carbon ratio) 19 , 
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or CH 4 - the latter will be used here. A missing preceding subscript for a pure species will be 
taken to mean that its El is in terms of its own mass, i.e. Y is identical to X. Clearly, Els can only 
be meaningfully summed if the subscripts match. 

It is worth emphasizing that the normalization used in defining the El does not account for 
propulsion efficiency and so does not represent the best overall reduction of net emissions from 
an engine. For instance, if a particular scheme reduced EI(no 2 )NO x by 10% but required 15% 
greater overall fuel flow rate to achieve the same thrust, there would be a penalty due to both poorer 
fuel consumption and greater emissions. The Environmental Protection Agency Parameter (EPAP) 
for landing-takeoff cycle conditions 29 represents a means of including the propulsion performance 
in the emissions estimation method. A similar approach could be used to compare engines in 
supersonic cruise operation for HSCT stratospheric impact assessment. 

Further, the total performance of the aircraft is not being considered - the most likely relationship 
being a greater aircraft weight resulting from emission reduction procedures 12-14 . Greater aircraft 
weight translates into more fuel burned and thus an increase in net emissions for a given El. Thus, 
the net emissions from an HSCT will depend on the total aircraft design, not simply on an El. 
On the other hand, combustor design will be optimized by lowering the important Els and, since 
improved combustor design will lower NO x emissions for any given propulsion system or airframe 
configuration, lowering Els is a key factor in improving the overall emissions performance of an 
aircraft. 

NO x ESTIMATION METHODS 

Extensive NO x emissions measurements have been performed over the years on many existing 
engines and engine components. Development work is proceeding towards a Tow NO x ’ HSCT 
combustor, but a complete prototype engine is still a long way from the test cells. Estimates 
for the expected emission performance of a given design are calculated by using semi-empirical 
correlations that have been developed and refined throughout the past emission measurement 
programs. These correlations account for the experimentally observed pressure and temperature 
dependences of the net NO x emission for a given combustor. Insofar as changes in combustor 
design do not alter the overall chemical kinetics, the correlations can be used to guide the combustor 
development efforts. 

The essentially square-root dependence of EI(no 2 )NO x on combustor pressure and the expo- 
nential dependence on combustor inlet temperature were used in correlations for CIAP work and 
current correlations have evolved from them. Over the range of combustor designs measured to 
date, these dependencies have been demonstrated to fit the observed data well, with a single mul- 
tiplicative coefficient reflecting a given combustor’s overall emission performance. For example, 
GE Aircraft Engines uses a correlation: 


EI(no 2 )NO x = 0.0986 


P 3 

1 atm 



n 

194.4 K 


Ho 


53. 2g H 2 0/kg dry air 


( 2 ) 


for their ‘dual annular’ low NO x combustor for the NASA ECCP 30 clean combustor program. 
In the correlation, p 3 is the combustor inlet pressure, and T 3 is the combustor inlet temperature. 
The effect of ambient humidity, H 0 in grams of water per kilogram dry air, is to reduce NO x 
production but at the low humidity levels encountered in the stratosphere, this correction will be 


5 



negligible. GEAE’s previous combustors used in the CF6 series engines are also approximated by 
this correlation with the multiplicative constant (0.0986) increased by 25% and a constant value 
of 2.2 added to the entire expression for the CF6-80C, and a 35% increase and 1.7 added for the 
CF6-50C. 

The coefficient of GEAE’s correlation is affected by the specific combustor design, including 
the residence time that a parcel of reacting fluid remains within the combustor. This particular effect 
is included parametrically in the correlations of Pratt and Whitney and NASA by referencing a 
representative combustor velocity. These two correlations are the same except that NASA’s does not 
correct for water vapor effects and both include a reciprocal reference velocity, V re f, dependence 
appropriate to the combustor flow and a dependence on T 4 , the combustor exit temperature, 
reflecting the temperature rise due to combustion heat release. 


EI(no 2 )NO x ~ 


pVt, ( t 3 


V ref 

EI(no 2 )NO x 


■exp 


H 0 


V r ef 


288 K 53. 2g H 2 0/kg dry air ) 
(NASA LeRC) 


(P&W) . 


p° 3 5 t 4 ( t 3 


■exp\ 


288 K 


(3) 

(4) 



Fig. GEAE’s correlation of EI(no 2 )NO x as a function of 10 atm < p 3 < 14 atm 
and 770 K < T 3 < 1000 K for the ECCP technology combustor. 
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The general functional dependence of these correlations is shown in the figure below, where 
EI(no,)NO x is plotted using the GEAE correlation for pressures from 10 to 14 atmospheres and 
temperatures between 770 K and 1000 K. Pratt and Whitney’s and NASA’s correlations have a 
similar functional form for reasonable values of T 4 , allowing for the combustion heat release and 
overall equivalence ratio. For T 4 ~ T 3 + 1100 K, a multiplicative constant can be chosen for 
the P&W/NASA correlation so that the temperature and pressure dependences agree with the GE 
correlation within 5-10%. 

For estimating EI (NO:! )NO x for the Lean, Premixed, Prevaporized (LPP) combustor, GEAE 
uses the correlation 31 : 

EI<no,)NO x = t„, exp( - 72.28 + 2.8 T^l - J , (5) 

where t rea is the combustor residence time in ms and T a di a is the adiabatic flame temperature in 
Kelvins. 

These correlations are used to estimate the emissions performance of a combustor for conditions 
that have yet to be tested. Indeed, the most promising candidate NO x reduction schemes have not 
yet been tested at the conditions expected in HSCT combustors. Promising experiments 32 indicate 
that advanced combustor concepts may achieve significant reductions at p 3 , T 3 appropriate to 
HSCT engine cycles. However this indication is based on experiments at less severe operating 
conditions extrapolated to the desired p 3 , T 3 . Future experiments will clarify the performance on 
this combustor approach at the anticipated HSCT conditions and further test the correlations for an 
as yet uncalibrated technology. 

EMISSION INDICES 

Emission indices were measured for several engines intended for supersonic flight during the 
1970’s 1 8,33 ,34 . The results of some measurements from this time period are tabulated below. 
Included are measurements of the one supersonic civilian transport that has operated commercially, 
the Olympus engine of the Concorde airplane. An important distinction must be made in comparing 
these emission indices to those calculated for the proposed HSCT engine. The increased propulsion 
efficiency inherent in the designs for the HSCT dictate higher combustor inlet temperatures (and 
pressures) than were used in the 1970’s engines. The NO x produced at these more severe conditions 
would be much higher than that measured in these earlier engines if the same combustor technology 
was used. Substantial improvements in the combustor design; improving the mixing, controlling 
the local equivalence ratio, and minimizing the residence time; are required to decrease the NO x 
emissions relative to these engines while simultaneously maintaining the higher operating pressure 
ratios, and thus combustor inlet temperatures, necessary for the required cycle efficiency. 

Also in the table included are emissions measurements for two afterburning military turbojet 
engines 18,33,34 . The measurements for the J58 engines were in preparation for the YF-12 in-flight 
NO x field measurements. The study of these engines demonstrated trends of Els versus altitude, 
fuel-air ratio, Mach number, and the presence or absence of afterburning for a specific engine. The 
optimized design of an engine designed to fly within a narrow range of flight conditions may not 
necessarily reflect all of these trends, however. 
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While it is significant that EI( N o 2 )NO x was reduced by a factor of ~ 2 during afterburning 
operation, this cannot be employed as a reduction method for HSCT propulsion. The reduction 
is due to the afterburning fuel being consumed^ at lower pressure antT temperature after expansion 
through the turbine, but this is coupled to increased El CO, El THC, and, most notably, substantial 
increases in specific fuel consumption and possibly more total NO x produced. Increases in spe- 
cific fuel consumption, in particular, could not be tolerated for sustained supersonic cruise in an 
economically viable commercial aircraft. 


TABLE 1 


Measured 1970’s Supersonic Engine Emission Indices 

Species 

Olympus 593 18 

YJ93-GE-3 18 

J58 33,34 

CO 

1.4 -4.6 

3.4-38.7 

1-120 

(no 2 )NO x 

17-20 

1.81 - 12.03 

2-20 

(No 2 )NO 

16-19 

1.83-10.69 


“ N0 2 

0.2 -1.3 



(ch 4 )THC 


0.03-0.8 

0.03 - 12 


Pratt and Whitney’s calculated emissions performance for their proposed Turbine Bypass 
Engine (TBE) 33 is included in Table 2. General emissions levels are given, but engine cycles have 
been optimized to different degrees for each case so, for example, the fact that EI( N q 2 )NO x is 
slightly greater for Mach 2.4 should not necessarily be construed as greater NO x emissions for that 
case. The changes in CO 2 and H 2 O emissions at Mach 3.2 relative to the lower speeds are due 
to a fuel change. The higher velocity flight conditions require a different fuel (JP-7 instead of the 
more commonly used Jet A) with a higher thermal stability and a concomitant increase in the H/C 
and H 2 O/CO 2 ratios. The soot/particulate El is calculated from the SAE Smoke Number, which is 

a quantitative measure of particulate emissions but is not directly based on the emitted particulate 

20 

mass . 
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TABLE 2 


Pratt and Whitney Estimates for En 

gine Emission Indices 35 

Species 

Mach 1.6 
18.3 km 

Mach 2.4 
19.8 km 

Mach 3.2 
21.3 km 

C0 2 

3165.5 

3165.5 

3117 

H 2 0 

1233 

1233 

1350 

CO 

1.1 

1.3 

1.9 

(NO 2 )N0 x 

5.3 

6.4 

5.1 

(no 2 )NO 

4.5 

5.4 

4.3 

no 2 

0.8 

1.0 

0.8 

so 2 

1.0 

1.0 

1.0 

soot/particulates 

0.02 

0.02 


(ch.)THC 

0.1 


0.2 

SAE Smoke Number 

18.4 

19.7 

14.0 

Fuel 

Jet A 

Jet A 

JP-7 


The calculated emissions performance of GE Aircraft Engines’ Variable Cycle Engine 36 (VCE) 
is quantified in Table 3. The NO x emissions were based on the correlation of eqn. 5. 


TABLE 3 


GE Aircraft Engines Estimates for Engine Emission Indices 3 ** 

Species 

Mach 2.4 

Mach 3.2 


16.8 km 

21.3 km 

C0 2 

3156 

3135 

h 2 o 

1240 

1290 

CO 

< 5 

<6 

(no 2 )NO x 

5.7 

7.0 

(NOj)NO 

5.1 

6.3 

no 2 


0.7 

so 2 

1.2 

1.2 

soot/particulates 

- 

- 

(ch 4 )THC 

<0.1 

<0.1 

SAE Smoke Number 

<1 

<1 

Fuel 

Jet A 

JP-7 


MTU Motoren- und Turbinen-Union in Munich Germany estimates emissions performance of 
HSCT engines 37 for both existing technology and expected future technology as indicated in Table 
4. 
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TABLE 4 


MTU Estimates for Engine Emission Indices 37 

Species 

Mach 2 
14.8 - 18.2 km 

Mach 3 
19.9 -23.6 km 

Advanced Proven Technology 

CO 

4.0 

4.0 

( N o 2 )NO x 

26-21 

19-15 

(CH 4 )THC 

0.4 

0.4 

Future Technology (15 to 20 years) 

CO 

2.6 

2.6 

(no 2 )NOx 

11-9 

8-7 

(ch 4 )THC 

0.3 

0.3 


In CIAP, consensus predictions were formulated 38 for exhaust emissions for anticipated emis- 
sion reduction technologies to be used in the SST program. Those predictions were for technological 
development that was never pursued to the demonstration of its potential. However, they serve as 
a point of reference for current predictions and are included in Table 5. The last column in Table 5 
gives the range of emission indices being considered for preliminary assessments. 


TABLE 5 


Estimates for HSCT Els in the Stratosphere 

Species 

CIAP 

Turbojet 38 

CIAP 

Duct-burning 38 

HSRP/AESA 
Current Range 

C0 2 



3120-3170 

•“ h 2 o 



1230- 1350 

CO 

3 

15-30 

1.1-6 

(no 2 )NOx 

3-14 

3-12 

5 - 50f 

(no 2 )NO 



4-45f 

S0 2 



0- 1 

soot/particulates 

0.02 

0.02 

— 

(ch 4 )THC 

0.1 -0.5 

3-10 

<0.1 -0.2 


f The upper value represents existing proven technology that is not likely to be used 
in its current form in a future HSCT aircraft. 


The range of C0 2 , and H 2 0 emissions are purely a function of fuel composition ^changes 
in these emissions are related to changes in fuel composition that increase its thermal stability as 
becomes necessary around Mach 3. S0 2 is also a function of the fuel itself - all the fuel sulfur 
is emitted as S0 2 . CO and THC emission estimates will have a broad range as long as both duct 
burning and unaugmented turbojet or low-bypass turbofan options are being considered, but are 
nonetheless bounded by fuel consumption constraints for an economically viable vehicle. 
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NO x emissions have the greatest range, as well as the greatest possibilities for control within 
the operational envelope that an HSCT will function. The upper limit is conservative, in the sense 
that it uses an estimate made with the GE correlation for the GEAE/NASA ECCP combustor 
for which measurements have been made. The low value of 5 is not a firm prediction of future 
technology, but rather an HSRP goal. These low values make use of extrapolations employing the 
correlations and basic laboratory data, sometimes using gaseous fuel experiments 31 , to estimate 
the best performance foreseeable at this time. It is hoped that, in the course of the current HSRP 
effort, these ranges in particular can be narrowed. 

AIRFRAME EMISSIONS 

The exhaust from the engines is the largest quantity of material left in the flight path by an 
airplane. However, the exhaust species that are directly involved in ozone chemistry are trace 
species. Because they comprise a small fraction of the total exhaust mass, the relative amounts of 
other minor emissions from the proposed HSCT need to be quantified to evaluate their potential 
contribution to the total aircraft emissions. These additional airframe emissions, as distinguished 
from the engine emissions contained in the exhaust, originate from 1) aircraft systems involved in 
controlling or supporting flight operations and 2) passenger systems providing services onboard. 
Estimates for these two classes of emissions for an HSCT will be made in the following subsections. 

The various fluids carried on the aircraft could potentially leak from the vehicle and find their 
way into the atmosphere. In making the following estimates, some worst-case assumptions will be 
made by assuming that an estimate for a fluid loss corresponds to an emission into the stratosphere. 
In fact, lost fluid could remain primarily within the aircraft, and that which escapes may do so when 
the plane is on the ground or during takeoff, landing, or while passing through the troposphere. The 
emissions escaping the vehicle at high altitude cruise could settle out of the stratosphere without 
vaporizing at the low temperatures (200-250 K) encountered. The estimates given are intended 
to provide a crude upper bound for the quantities of these species that might be released in the 
stratosphere. If any of these estimates suggest that an airframe emission species is important to 
stratospheric chemistry, more refined quantitation would be required to model its effects. 

Since the exhaust emissions are reported as emission indices (Els) in grams of emitted species 
per kilogram fuel burned, a comparison of those exhaust Els and the aircraft emissions requires a 
correspondence relating these two quantities. Aircraft emissions are not caused by the consumption 
of fuel in any direct sense, but rather are associated with 1) routine operation and thus flight time 
or number of flights, or 2) random, unusual events with some small frequency of occurrence. 
Conversion of these aircraft emissions to an equivalent emission index (EEI) will be made by 
using an estimate of the fuel bum rate appropriate for the generalized design of a Mach 2.4 HSCT 
aircraft and making the conversion from time to fuel mass. 

Aircraft Systems 

All the mechanical and electrical systems needed to control an HSCT in flight have not 
yet been designed in detail, so the amounts of mat erials carried onboard or used during operation 
cannot be specified. Corresponding quantities from existing aircraft cannot be substituted without 
qualification, since particular systems can be very configuration dependent; beyond that, the HSCT 
is a substantially different vehicle from a subsonic aircraft, thus its systems may differ significantly 
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from those aircraft for which quantities are available. On the other hand, the types of systems 
needed are basically the same as those currently in use, and available information on existing 
aircraft can be used to provide crude estimates of emissions due to these sources. 

The most ample fluid carried is, of course, the engine fuel. The economical use of fuel is a 
driving criterion in designing the aircraft, so little unbumed fuel is allowed to escape the aircraft. 
Fuel is dumped very occasionally if too much is being carried for a safe landing, but this is an 
emergency measure carried out when approaching a landing and the aircraft is below about 2 km. 
Thus, fuel is not dumped in the stratosphere. Fuel tanks may be vented during flight, releasing 
some vaporized fuel, and this may be a bigger factor for a supersonic aircraft with significant skin 
heating than it is for subsonic aircraft, but no estimate for such vapor loss has yet been obtained. 

Hydraulic fluid is used to control mechanical devices in the aircraft, and regular maintenance 
requires replenishing the hydraulic system reservoirs. Boeing 39 estimates that regular seepage or 
leakage amounts to 20 gallons (70 kg) of hydraulic fluid (composed of phosphate esters) lost from 
the hydraulic systems per year per airplane. In addition, aircraft are designed with redundant 
hydraulic systems to ensure control should one system fail. Failures are rare events, even so; an 
estimate of an average loss rate by Boeing 39 is 1 gallon (3 kg) per year per airplane due to system 
failure. (This is based on a loss of 6 gallons of fluid on average in arf event that occurs 50 times in 
10 6 flight hours, as extracted from reliability reports for existing subsonic commercial aircraft, and 
33 1 1 flight hours per year per airplane for these subsonic transports.) While it will be assumed that 
these 70 kg per year per airplane (or 2 x 10 4 kg/ 10 6 flight hours) of phosphate esters are deposited 
in the stratosphere, for the low vapor pressure hydraulic fluids especially, it is not likely that this 
emission could appear in the stratosphere as a chemically reactive species. 

Oil is used as a lubricant for the engines, as well as some auxiliary devices; infrequently 
some oil is discharged due to accidental leaks. Some oil/air mist may be vented from the oil 
reservoir tanks 40 . Some oil is consumed within the engine, as well - some of it being burned in the 
combustor, some exiting with the exhaust. The emissions due to oil consumed within the engine 
would be quantified in a measurement of the total exhaust emissions and thus really represent 
an engine emission in d iscernible from the emissions due to fuel consumption. It is worth noting 
however, that as an engine ages, changes in oil consumption may result in changes in exhaust 
emissions as well. 

Accidental leakage of oil, like hydraulic fluid, will not necessarily leave the aircraft, nor will 
that which leaves necessarily remain in the stratosphere. However, as an upper bound, Boeing 39 
quotes a value of 5 incidences in 10 6 flight hours releasing on average 20 kg of oil. This amounts 
to a loss of 0.3 kg per year per aircraft of a fluid composed primarily of long chain aliphatic 
hydrocarbons with trace additives of phosphate esters and metal compounds. 

Other less quantifiable emissions 41 could arise from the venting of the lead-acid and/or 
nickel-cadmium batteries used for reserve power; this would result in negligible amounts of H 2 and 
O 2 being released. (Emissions due to mists of KOH and H 2 S0 4 electrolytes would be even smaller 
and not likely to leave the aircraft.) Auxiliary power units (APUs) are small gas-turbine engines 
that would not be used in cruise typically and would emit much smaller amounts of exhaust not 
unlike the main engines. Windshield washer fluid is more likely to be used at lower altitudes, and 
the glycol solution used is unlikely to amount to even a small fraction of the ‘blue ice’ emissions 
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(see below). Finally, novel anti-icing techniques may make use of glycol solutions that would be 
carried on board, although not used in the stratosphere above hazardous clouds, but emissions due 
to this fluid, too, are not expected to amount to any more than the ‘blue ice’ fluids at stratospheric 
altitudes. 

Passenger Systems 

A variety of systems within the airplane’s cabin are devoted to serving the passe’ „ers on 
board. Several refrigeration systems are typically available for air conditioning, refrir ition of 
food and drink in the galley, and sometimes cooling of avionics (more correctly an aircr ft system 
need). Until replacement refrigerants are found, these systems will contain CFCs, and any leakage 
that exits the aircraft in the stratosphere will be depositing CFCs where they are photochemically 
active. Douglas 41 estimates that 0.2 kg per year per airplane is typically used to maintain a 
representative aircraft’s refrigerant systems. 

Galley cooling and avionics cooling do not necessarily require CFC-based refrigeration 
(recirculating liquids are currently employed in some situations), so the amount of refrigerant 
aboard a yet-to-be-developed airplane could be minimized, if direct deposition in the stratosphere 
must be minimized. In a similar vein, fire extinguishers on current aircraft contain halons 41 - a 
source of bromine and chlorine when photochemically activated in the stratosphere. Whether there 
is any greater risk due to stratospheric flight with CFCs and halons (or their replacements) than 
flight in the upper troposphere with existing aircraft would determine the need for replacement 
or closer control of leaks and accidental release on board an HSCT. For long-lived halocarbons, 
release in the lower stratosphere is not substatially different from release at the surface. 

Lavatory holding tanks have received some media attention over the past several years as a 
result of mishaps attributed to ‘ blue ice’ - lavatory fluid that has leaked from the holding tank due to 
faulty maintenance, frozen in the cold atmosphere, and then broken free. This has caused damage 
to the aircraft by being ingested in the engine or to property on the ground from impact. From the 
point of view of safety, this is a major problem and efforts have been directed at totally eradicating 
any possibility of leakage. As far as emissions are concerned, the low (and probably decreasing) 
frequency of occurrence and the likely substantial descent before complete vaporization suggest 
that lavatory fluids will not contribute to the total aircraft emissions in the stratosphere. Based on 

1 incident involving 30 kg per 10 6 flight hours, Boeing 39 estimates that less than 1 cup of fluid is 
lost per year per aircraft using current reliability data. 

For subsonic aircraft, oxygen service is available for passengers ‘in the event of sudden 
loss of cabin pressure’. While this service is seldom used, some system leakage might occur 
from gaseous systems 42 (used for crew and perhaps passengers for high altitude flight). Chemical 

0 2 generation packs 43 as used in many subsonic airplanes would only emit during use. Purity 
requirements for passenger respiration units limit Cl, CO, and CO 2 levels to 0.1, 20, 1000 ppmv 
respectively, so even if used, only oxygen with trace impurities is released. However oxygen is a 
major atmospheric constituent and even a major release of gas from the emergency oxygen system 
could not affect stratospheric chemistry. The substantially lower static pressures in the stratosphere 
necessitate more sophisticated depressurization safety measures for HSCTs. Potential emissions 
from such safety systems must be evaluated, particularly if they use fluids not already used in other 
systems. 
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Calculation of Equivalent Emission Indices for Airframe Emissions 

To compare the amounts of these airframe emissions with those from the engine, an equiv- 
alent emission index will be defined. Simply, the average emission amount in grams released in 
a given flight time will be divided by the total amount of fuel used during cruise in the same time 
period, i.e. 

average emission rate of X g/s 
fuel burn rate kg/s 

Using the total of 33 1 1 flight hours per year per aircraft for the existing subsonic aircraft (as used 
above) to arrive at a mean emission rate from the subsonic reliability data and taking a fuel bum rate 
of 9 kg/s as repr esentati ve for Mach 2.4 flight, th e emission estimates noted above can be converted 
to EEIs for HSCT cruise and are listed in Table 6. Note that a similar value of 3100 stratospheric 
flight hours is considered representative for HSCT flight 39 as well. Again, it is assumed that all the 
lost fluids are emitted and vaporized in the stratosphere, which in many of these cases is not likely 
to happen and so these quantities should be interpreted as upper bound estimates. 


TABLE G 


Upper Bound Estimates for Airframe Emissions in the Stratosphere 

Emission 

kg / 10 b flight hours 

Equivalent Emission Index 

Aircraft Systems 


lubrication oil burst 

1 xlO 2 

3 x 10 -6 

hydraulic fluid burst 

' 1 xlO 3 

3 xlO" 5 

hydraulic fluid leaks 

2 xlO 4 

-r 

1 

O 

1 4 

X 

fuel dump 

n/af 

n/af 

fuel vent 

no estimate 

no estimate 

Passenger Systems 


refrigerants 

1 xlO 2 

4 xlO -6 

‘blue ice’ 

30 

1 xlO -6 


f Not Applicable. Fuel is not dumped at stratospheric altitudes. Corresponding 
tropospheric quantities would be 3 x 10 6 kg / 10 6 flight hours, EEI = 0.1 39 . 

Referring to Tables 2-4, it is apparent that the largest stratospheric airframe emission is 
substantially smaller than the smallest engine emission and several orders of magnitude smaller 
than NO y emissions. Recalling that the airframe emissions are upper bound estimates, it appears that 
these airframe emissions can be neglected relative to those from the engine. One notable exception 
might be refrigerants and gaseous fire extinguishers; if these species are highly photochemically 
active, the effects of direct deposition in the stratosphere need to be assessed. If there are no highly 
reactive catalytic effects of any of the other airframe emissions, it appears that their contribution 
to the net aircraft emissions may be neglected. 
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CONCLUDING REMARKS 

Estimates have been provided for the emissions from a proposed HSCT in stratospheric flight. 
The emissions from the aircraft are expected to be dominated by the effluents from the engine and 
the emissions from the airframe are probably negligible. Current estimates for the NO x emissions 
have been calculated by engine companies for the proposed engine cycles using correlations that 
have been developed over the past several decades of NO x reduction programs. These estimates 
will have to be validated with future tests of components and full engines as development proceeds. 

More complete measurements of particulate size distributions and number densities and of 
aerosol-active species will be required if heterogenous chemistry is to be accurately modeled, 
both in the wake of individual airplanes and globally. New measurements of NO y speciation will 
be necessary to completely catalog the ‘odd- nitrogen’ emissions and to resolve uncertainties in 
measuring high temperature, high velocity exhaust gases with a sampling probe. 
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Chapter 2 

Plume and Vortex Regime Modeling of Mixing and Chemistry* 

SUMMARY 

The environmental perturbations caused by the exhaust of a High Speed Civil Transport (HSCT) 
depend on the deposition altitude and the amount and composition of the emissions. The chemical 
evolution and the mixing and vortical motion of the exhaust need to be analyzed to track the exhaust 
and its speciation as the emissions are mixed to atmospheric scales. Elements of an analytic model 
of the wake dynamical processes are being developed to account for the roll-up of the trailing 
vorticity, its breakup due to the Crow instability, and the subsequent evolution and motion of 
the reconnected vorticity. The concentrated vorticity is observed to wrap up the buoyant exhaust 
and suppress its continued mixing and dilution. The chemical kinetics of the important pollutant 
species will be followed throughout the plume and wake. Initial plume mixing and chemistry 
are calculated using an existing plume model (SPF — Standard Plume Flowfield) with additional 
H/C/O, 0H/S0 2 , and NO y chemistry and equilibrium H 2 0 condensation included. The species 
tracked include those that could be heterogeneously reactive on the surfaces of the condensed solid 
water (ice) particles when condensation occurs and those capable of reacting with exhaust soot 
particle surfaces to form active contrail and/or cloud condensation nuclei (ccn). 

INTRODUCTION 

The Anglo-French Concorde has shown that commercial supersonic air transportation is fea- 
sible. However, since it began commercial operations in the mid-1970s, the Concorde has served 
a low-volume, premium-fare market with its dozen-plus fleet. Technological advances in the in- 
tervening years have created the potential for a large fleet of more economical High Speed Civil 
Transports (HSCTs) that could serve a large commercial market, 1,2 including trans-Atlantic and 
Pacific Rim routes in particular. With this possible development, attention has refocused 3 on the 
environmental impact of a much larger fleet of airplanes operating in a flight regime close to that 
of the Concorde. 

The potential effects of these aircraft on the chemical balance in the stratosphere — and on 
stratospheric ozone levels in particular — are a primary environmental concern. Atmospheric 
modelers 4,5 have already begun to make estimates of the effects of the exhaust effluents of HSCTs 
on stratospheric chemistry. In doing so, a range of estimated emissions from the as-yet-undeveloped 
HSCT engines are deposited directly at representative altitudes and latitudes in their models. It is 
the purpose of the present study to provide initial estimates of chemical and physical changes that 
may occur between the exhaust leaving the engine and its deposition as a quiescent air mass in the 
stratosphere. 

The exhaust from the large engines propelling a 250 + passenger supersonic airplane will 
be introduced into the ambient surroundings as a supersonic coflowing jet. Recent studies of 
supersonic mixing 6,7 indicate that supersonic convective Mach numbers suppress mixing in free 
shear layers, an effect that will need to be accounted for in following the initial exhaust dilution. 
As this initial dilution is occurring, the vortex sheet shed by the lifting surfaces of the airplane is 

* This chapter has been presented as AIAA paper 91-3158 at the Aircraft Design and Systems Meeting, 23-25 
September 1991, Baltimore, MD. Also accepted for publication in J. Aircraft . 
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rolling up and being concentrated in a trailing vortex pair by mutual induction of the distributed 
vorticity. This vorticity field will begin to affect the diluting exhaust stream and will subject the 
exhaust to local conditions that are sufficiently different from ambient that significant changes in 
stratospheric exhaust deposition could result. The vorticity field behind a highly loaded supersonic 
aircraft configuration can differ dramatically from that of typical subsonic transports. 

In the Climatic Impact Assessment Program of the early 1970s, a related study 8 was performed 
using contrail data collected under that program and then-current understanding of chemistry, 
condensation physics, and vortex wake development. Just as supersonic flight technology has 
made great bounds since that time, revolutionary advances have recently been made in atmospheric 
chemistry. Great strides in the understanding of the stability and development of vortex wakes 
were being made concurrently with the CIAP program and continued since then and were not fully 
incorporated into that early study. Thus, while not starting completely anew, the present study 
represents a significant departure from the CIAP study and will provide a more comprehensive 
assessment of the chemistry and fluid physics occurring in the vortex wake. 

The elements of an analytic model of these processes must include the dynamics of the vortex 
wake, from roll-up to breakup and dissipation; its interaction with the diluting hot exhaust; the gas- 
phase exhaust species chemistry; the heterogeneous chemistry occurring on condensed aerosols 
and soot particles; and the condensation physics governing the creation, or lack thereof, of the 
condensed aerosols needed for heterogeneous reactions. The stratospheric impact of the net result of 
all these elements will be felt through their displacement of the exhaust through vortical or buoyant 
transport 9 or possible sedimentary removal of reactive species, as in polar stratospheric cloud (PSC) 
denitrification (and dehydration) of the polar stratosphere. 10 ’ 11 Any major chemical repartitioning 
of the exhaust species that is not washed out by subsequent photochemical conversions may 
significantly affect the composition of the chemical inputs to the global models. An understanding 
of the structure, motions and persistence of the vortex wake will also be necessary for any attempt 
to measure the exhaust from an HSCT in flight and for evaluating the wake hazard created for 
following aircraft. 

The analysis reported below represents the development of these several elements of a compre- 
hensive exhaust plume/wake model for supersonic aircraft in cruise conditions. Gas-phase nozzle 
chemistry and plume chemistry are described in the next section and have been applied to under- 
standing the engine exit plane conditions and the subsequent chemical evolution and condensation 
physics in the initial exhaust mixing region using extensions of a standard aircraft “plume" code. 
This chemistry will be applied in the vortex wake regime as well, as a refined model of its dynamics 
is finalized. The analytical basis for the vortex wake model is presented in the following section. 
The integration of these elements and requirements for further development are discussed in the 
last summary section. 
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TABLE 1. 

Reaction List for HSCT Internal Engine Flow Model 

reaction rate constant (cm 3 molecule - 1 s — *) 

k = A e (~ E / RT ) /T n 

A N E/R (K) 


2H + M = H 2 + M 
20 + M = 0 2 + M 
H + O + M = 0H + M 
OH + H + M = H 2 0 + M 
CO + O + M = C0 2 + M 
20H = H 2 0 + O 
OH + H 2 = H 2 0 + H 
OH + H = H 2 + O 
OH + O = H + 0 2 
OH + CO = C0 2 + H 
H + 0 2 + M = H0 2 + M 
H + H0 2 = 20H 

ho 2 + h = h 2 + o 2 

H0 2 + H 2 = H 2 0 2 + H 
2 OH + M = H 2 0 2 + M 
2H0 2 = h 2 o 2 + o 2 
h 2 o 2 + oh = h 2 o + ho 2 

NO + o = 0 2 + N 
O + N 2 = NO + N 
NO + H = OH + N 
N0 2 + H = NO + OH 
N0 2 + O = NO + 0 2 
n 2 o + O = 2NO 
N 2 0 + H s N 2 + OH 

no 2 + m = no + o + m 
2N0 2 = 2N0 + 0 2 
H0 2 + OH = H 2 0 + 0 2 
H0 2 + O = OH + 0 2 
H0 2 + M = H + 0 2 + M 
CO + 0 2 = C0 2 + O 
CO + H0 2 = C0 2 + OH 
H 2 0 2 + H = H 2 0 + OH 
h 2 o 2 + O = OH + ho 2 
co 2 + o = CO + o 2 

N0 2 + OH + M = HNO3 + M 


1 . 50 X 10 29 

- 1.3 

0.0 

5 . 20 X 10 -35 

0.0 

- 900.0 

— 29 

1 . 30 X 10 

- 1.0 

0.0 

6 . 10 X 10 -26 

- 2.0 

0.0 

1.70 X 10 — 33 

0.0 

1510.0 

3 . 50 X 10 -16 

1.4 

- 200.0 

1 .06 X 10 — 17 

2.0 

1490.0 

8 . 10 X 10 -21 

2.8 

1950.0 

7 . 50 X 10 -10 

- 0.5 

30.0 

2 . 80 X 10 -17 

1.3 

- 330.0 

1 . 77 X 10” 30 

- 1.0 

0.0 

2.80 X 10 — 10 

0.0 

440.0 

1 . 00 X 10 -10 

0.0 

1070.0 

5 . 00 X 10 -11 

0.0 

13100.0 

1.60 X 10 — 33 

- 3.0 

0.0 

3 . 00 X 10 -12 

0.0 

0.0 

2 . 90 X 10 -12 

0.0 

160.0 

6 . 31 X 10 -15 

1.0 

20820.0 

3.02 X 10 — 10 

0.0 

38370.0 

4 . 37 X 10 -10 

0.0 

25370.0 

5 . 76 X 10 -10 

0.0 

740.0 

1.66 X 10 — 1 1 

0.0 

300.0 

1.15 X 10 — 10 

0.0 

13400.0 

1.26 X 10 — 10 

0.0 

760.0 

1.82 X 10 — 08 

0.0 

33000.0 

3.31 X 10 -12 

0.0 

13500.0 

2.40 X 10 — 08 

- 1.0 

0.0 

2.90 X 10 — 1 1 

0.0 

- 200.0 

2 . 00 X 10 -05 

- 1.2 

24363.0 

4.20 X 10 — 12 

0.0 

24000.0 

2.50 X 10 — 10 

0.0 

11900.0 

4 . 00 X 10 -11 

0.0 

2000.0 

1.60 X 10 — 17 

2.0 

2000.0 

2 . 80 X 10 -11 

0.0 

26500.0 

— 22 

2 . 20 X 10 

- 3.2 

0.0 



PLUME MIXING AND CHEMISTRY 

Internal Engine Flow/Finite- Rate Kinetics Model 

The internal flow for a hypothetical Mach 2.4 engine design was modeled under realistic engine 
operating conditions using the ODK PACKAGE code 12 for reacting flows. Emphasis was given to 
estimating nonequilibrium concentrations of OH and HNO3 at the nozzle exit plane since emission 
levels for these species are not currently being measured by propulsion laboratories and they 
were expected to strongly influence subsequent heterogeneous and homogeneous exhaust plume 
chemical kinetics. 

The finite-rate kinetics for the internal engine flow were simulated using the reaction set 
"consisting of 35 reversible reactions given in Table 1. These reactions were selected from kinetic 
data bases for hydrocarbon and nitrogen combustion chemistry. The hydrocarbon combustion 
reaction rate parameters were taken from a recent compilation 13 of kinetic rate data for methane 
combustion. Rate parameters for the reactions describing the nitrogen combustion chemistry were 
taken from the recent review by Miller and Bowman. 14 



TABLE 2. 



Internal Engine Flow Properties* 


combustor 

turbine turbine mixing chamber 

exit plane 

equil. 

equil. ODK equil. ODK 

equil. ODK 


p (atm) 

13.609 

11.3 

11.3 

1.18 

1.20 

0.042 

0.042 

T (K) 

2050.0 

1398.0 

1383.0 

1046.0 

1049.0 

430.0 

434.0 

time (s) 

NA 

0.0033 

0.0033 

0.006 

0.006 

0.0009 

0.0009 

Mach no. 

NA 

2.17 

2,17 

0.26 

0.26 

2.85 

2.84 


CO 

1.20 

X 

10" 

-4 

1.52 

X 

10 

-7 

2.44 

X 

10 

-5 

2.02X10" 

-11 

1.07 

X 

10" 

-5 

5.67X10“ 

-31 

1.09 

X 

10" 

-5 

C0 2 

6.84 

X 

10" 

-2 

5.94 

X 

10 

—2 

5.93 

X 

10 

—2 

3.22 X 10 

-2 

3.24 

X 

10" 

-2 

3.22X10’ 

-2 

3.24 

X 

10" 

-2 

H 

3.26 

X 

10" 

-6 

1.21 

X 

10 

-9 

1.54 

X 

10 

-6 

3.52X10" 

-14 

1.49 

X 

10" 

-9 

0.0 


6.12 

X 

10" 

-9 

HNO 3 

3.28 

X 

10" 

-10 

2.09 

X 

10 " 

-11 

6.24: 

x: 

L0" 

-10 

3.06X10" 

-11 

4.46 

X 

10" 

-9 

4.78X10" 

-12 

4.07 

X 

10“ 

-8 

ho 2 

1.61 

X 

10" 

-6 

1.93 

X 

10 

-8 

6.97 

X 

10 

-7 

3.05X10" 

-10 

1.04 

X 

10 " 

-7 

1.03X10" 

-20 

9.77 

X 

10“ 

-8 

h 2 

2.84 

X 

10" 

-5 

7.91 

X 

10 

-8 

7.87 

X 

10 

-6 

2.70X10" 

-11 

1.07 

X 

10“ 

-6 

4.17X10" 

-28 

1.03 

X 

10“ 

-6 

h 2 o 

7.61 

X 

10" 

-2 

6.62 

X 

10 

-2 

6.61 

X 

10 

-2 

3.57X10' 

-2 

3.59 

X 

10* 

-2 

3.57X10' 

-2 

3.59 

X 

10“ 

-2 

h 2 o 2 

1.90 

X 

10 " 

-7 

2.61 

X 

10 

-9 

1.65 

X 

10 

-8 

7.23X10" 

-11 

2.56 

X 

10“ 

-8 

3.66X10" 

-19 

2.44 

X 

10“ 

-8 

NO 

5.99 

X 

10 " 

-3 

5.54 

X 

10 

-4 

5.20 

X 

10 

-3 

4.85X10' 

-5 

2,81 

X 

10" 

-3 

1.66X10" 

-11 

2.81 

X 

10“ 

-3 

no 2 

2.17 

X 

10" 

-5 

3.16 

X 

10 

-6 

7.29 

X 

10 

-6 

1.76X10' 

-6 

8.72 

X 

10“ 

-6 

1.70X10' 

-9 

8.73 

X 

10“ 

-6 

n 2 

7.47 

X 

10" 

-1 

7.54 

X 

10 

-1 

7.52 

X 

10 

-1 

7.66X10' 

-1 

7.68 

X 

10“ 

-1 

7.66X10' 

-1 

7.68 

X 

10“ 

-1 

n 2 o 

1.08 

X 

10" 

-6 

3.31 

X 

10 

-8 

7.89 

X 

10 

-7 

3.45X10" 

-9 

4.30 X 

10* 

-7 

8.74X10" 

-16 

4.30 

X 

10" 

-7 

0 

7.95 

X 

10 " 

-5 

2.76 

X 

10 

-7 

3.23 

X 

10 

-5 

2.17X10" 

-10 

1.56 X 

10’ 

-7 

1.22X10" 

-27 

1.55 

X 

10" 

-7 

OH 

9.41 

X 

10' 

-4 

1.93 

X 

10 

—5 

2.09 

X 

10 

— 4 

1.39X10" 

-7 

4.21 

X 

10' 

-6 

8.13X10" 

-19 

3.40 

X 

10' 

-6 

o 2 

9.23 

X 

10" 

-2 

1.11 

X 

10 

-1 

1.08 

X 

10 

-1 

1.56X10" 

-1 

1.56 

X 

10' 

-1 

1.57X10' 

-1 

1.56 

X 

10" 

-1 


* species concentrations are in mole fraction 


To simulate the internal engine flow, equilibrium calculations were first performed to determine 
the chemical composition of the flow from the combustor as a function of temperature and pressure. 
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ODK calculations were then performed in three stages to characterize the chemical and physical 
properties of the flow through the turbine, mixing chamber and nozzle. For each stage, the 
pressure and temperature at the initial and end axial points were constrained to be consistent with 
representative GE engine parameters. 15 In addition, when necessary, the mixture was diluted to 
account for the mixing in of bypass air. Nozzle exit plane species mole fractions obtained by this 
method are compared against equilibrium predictions in Table 2. 

Although the one dimensional exit plume results shown in Table 2 can only approximate the 
chemistry and flow dynamics of anticipated HSCT engines, they do confirm that significantly non- 
equilibrium levels of CO, NO, N0 2 , HN0 3 and OH can be anticipated at the engine nozzle exit 
plane. While the engine development community has long accepted the non-equilibrium nature of 
CO and NO x emissions, there has been less appreciation of the high (> ppm) levels of OH emitted 
by gas turbine systems. The predictions presented here are consistent with the optical absorption 
measurements performed on the YJ93-GE-3 under stratospheric cruise conditions by McGregor et 
a i. 16 during the CIAP program. In that study, mole fractions of exit plane OH exceeded 2 x 10 -5 
and 4 x 10 -5 for cruise power settings equivalent to Mach 2.6 operation at 66 kft and Mach 2.0 at 
55 kft. 

HSCT Plume Flow/Finite-Rate Chemical Kinetics Model 

The standard plume flowfield code, SPF-2, 17 was used to characterize the chemical and physical 
properties of the combustion exhaust plume associated with an advanced high altitude supersonic 
aircraft. Model simulations used GE parameters 15 for a hypothetical Mach 2.4 engine to specify 
initial exhaust emission parameters and 1990 atmospheric models 18 to specify the ambient atmo- 
spheric parameters. The SPF code was modified to include an initial estimate of the degree of 
H 2 0 condensation based on thermodynamic equilibrium. Initial examination of the plume chem- 
istry focused on the evolution of NO x , including an evaluation of the impact of ozone oxidation 
mechanisms and the percentage conversion to HN0 3 , as well as SO x oxidation cycles. 

Initial Conditions The species concentrations and thermodynamic properties used as 
initial conditions in model simulations of the plume chemistry and mixing are given in Table 
3. This includes ambient free stream conditions and the initial conditions for the exhaust at the 
nozzle exit plane. Exit plane species distributions are based either on reported emission indices 
(g equiv/kg fuel) or on ODK calculations for the internal engine flow. An exit plane OH mole 
fraction of 1 x 10 -5 representing a rough mean between calculated ODK value shown in Table 2 
and the measurements of McGregor et ai. 16 was adopted for these plume calculations. For NO x , 
an emission index 19 EI(no 2 )NO x of 5 was used 20 assuming 10% N0 2 and 90% NO on a molar 
basis. Ambient free-stream conditions are given for three representative cases. These correspond 
to the atmospheric conditions appropriate for January 1 and June 1 at a latitude of N47 and June 
1 at N85. In all three cases, the altitude was ~18 km. Ambient conditions were obtained from 
AER’s 1990 atmospheric chemistry data base. 18 
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TABLE 3. 

Initial Conditions for HSCT Plume Flowfield Model* 



exhaust 


free stream 




case 1 

case 2 

case 3 

latitude 


N47 

N47 

N85 

altitude(km) 


18.4 

17.85 

18.3 

date 


June 15 

January 1 

January 1 

p(atm) 


0.0734 

0.0712 

0.0573 

T (K) 

5.61 X10 2 

2.19 X 10 2 

2.17 X 10 2 

2.05 X 10 2 

V (ft/s) 

4.30X10 3 

2.34 X 10 3 

2.32 X 10 3 

2.26 X 10 3 

CO 

2.37X10 -5 

1.99X10 -8 

1.98 X 10 ~ 8 

1.69 X 10 — 8 

co 2 

3.17 X 10 — 2 

3.50 X10 -4 

3.50X10 -4 

3.50 X 10 — 4 

H 

1.00X10 -7 

1.43X10 -20 

2.81 X10 -21 

0.0 

h 2 

0.0 

5.18X10 -7 

5.19 X 10 — 7 

5.20 X 10 _ 7 

h 2 o 

3.02X10 " 2 

4.20X10 -6 

4.20 X10 -6 

4.86X10 ~ 6 

H0 2 

0.0 

1.02X10 -12 

3.30X10 -13 

0.0 

h 2 o 2 

0.0 

3.51 X 10 " 12 

2 . 11 X 10 -12 

0.0 

n 2 

0.779 

0.790 

0.790 

0.790 

NO 

4.32X10~ 5 

2.43X10 -10 

5.87X10 -11 

0.0 

no 2 

4.80X10“ 6 

8.50 X 10~ 10 

4.29X10 -10 

0.0 

O 

0.0 

6.06X10 -13 

3.11X10 -13 

0.0 

OH 

1.00X10 -5 

1 .27 X 10 — 13 

2.76 X 10 — 14 

0.0 

0 2 

0.159 

0.210 

0.210 

0.210 

O 3 

0.0 

2.48 X 10~ 6 

2.46 X10 -6 

3.02X10 ~ 6 

HNO 3 

0.0 

3.71 X10 -9 

3.32 X 10 ~ 9 

7.39 X 10 — 9 

NO 3 

0.0 

1.23X10 -12 

1.71 X 10 — 12 

0.0 

n 2 o 5 

0.0 

3.57X10 -10 

6.09X10 -10 

0.0 

S0 2 

6.91 X 10 -6 

0.0 

0.0 

0.0 


* species concentrations are in mole fraction 


. : , , TABLE 4. _ 

SPF Equilibrium H 2 O Condensation Predictions 



case 1 

case 2 

case 3 

latitude 

N47 

N47 

N85 

altitude(km) 

18.4 

17.85 

18.3 

date 

June 15 

January 1 

January 1 

T(K) 

219.3 

216.7 

205.2 

p(atm) 

0.0734 

0.0712 

0.0573 

SPF Equilibrium H 2 0 




Condensation predictions 

no 

no 

yes 

Condensation algorithm 




predictions (Ref. 21 ) 

never 

uncertain 

always 



H 2 0 Condensation The SPF-2 code was modified to predict the degree of H 2 0 con- 
densation during the plume expansion as part of an initial investigation into the potential role of 
heterogeneous kinetic processes in affecting plume chemistry. Condensation was treated using 
the standard vapor pressure test for thermodynamic equilibrium. It was assumed that the any 
condensed H 2 0 particulates that formed would not significantly modify the flow mixing dynamics. 
A summary of the SPF equilibrium H 2 0 condensation predictions are given in Table 4 for the three 
case runs described by the free stream ambient conditions in Table 3. For comparison, condensation 
predictions based on Appleman’s equilibrium, fully-mixed condensation algorithm 21 are also listed 
in Table 4. In general, the agreement between these two approaches is good. 

Plume Chemistry 

The reaction set used to simulate the exhaust plume finite-rate chemical kinetics is listed in 
Table 5. Rate parameters for these reactions are based on a 1990 NASA evaluation 22 of the 
chemical kinetics data for use in stratospheric modeling. The reaction set includes reactions to 
describe NO x /C >3 reactions and S0 2 oxidation kinetics. 

Representative model results are presented in Table 6 which lists the plume mole fractions 
for nitrogen and sulfur oxides, nitric acid, hydroxyl radical and ozone on the plume axis at axial 
points 500, 1000 and 2000 feet downstream of the nozzle exit plane. The speciation given Table 
6 is based on model runs using ambient atmospheric conditions appropriate for June 15 at a N47 
latitude (case 1 in Table 3). However, similar results were obtained for model runs with ambient 
conditions corresponding to January 1 at the same latitude (case 2 in Table 3). The only difference 
was a 1.7 increase in the mole fraction of N 2 Os. 

Mole fraction ratios for N0 2 /NO, HN0 3 /NO x and S0 3 /S0 2 are also given in Table 6 where, 
in this case, NO x corresponds to the sum of NO and N0 2 mole fractions. Based on this data, 
the SPF model calculations indicate that approximately 11% of the initial S0 2 is oxidized to 
S0 3 . NO is oxidized to N0 2 by entrained ambient 0 3 . N0 2 , in turn, is oxidized primarily by 
exhaust OH. The N0 2 /NO ratio is maintained at 0.10-0.11 while approximately 5% of the NO x 
is converted to HN0 3 . Previous studies by Hoshizaki et a l. s concluded that about 10-20% of the 
NO x emitted is transformed to HN0 3 in the jet plume region, i.e. 1-10 seconds after the exhaust 
is emitted. Although this represents a significantly larger percentage conversion than predicted in 
SPF-2 simulations, the discrepancy can be explained in terms of differences in the rate parameters 
used. The dominant mechanism for formation of HN0 3 is 

Iq-22 

N0 2 + OH + M = HN0 3 +M, k=2.2x-^ T . 

In the early work of Hoshizaki et al. 8 this reaction was treated as a simple bimolecular reaction 
with a rate constant (k = 4.8 x 10 -12 ) orders of magnitude larger than that adopted in the 1990 
NASA evaluation 22 of the best fit to experimental data over the temperature range of interest. 
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TABLE 5. 

Reactions List for Plume Chemistry/Mixing Model 


reaction rate constant (cm 3 molecule *s *) 

k = Ae(- E ' RT ^/T N 



A 

N 

E/R (K) 

H + 0 2 + M = H0 2 + M 

5.2X10 -28 

1.6 

0.0 

OH + O = H + 0 2 

2.2X10 -11 

0.0 

240.0 

OH 4- H 2 = H 2 0 + H 

S.5X10 -12 

0.0 

-4000.0 

OH + OH = H 2 0 + O 

4.2 X 10~ 12 

0.0 

-480.0 

OH + OH + M = H 2 0 2 + M 

_?Q 

6.6X10 ™ 

0.8 

0.0 

OH + H0 2 = H 2 0 + 0 2 

4 . 8 X 10 -11 

0.0 

500.0 

OH + H 2 0 2 = H 2 0 + H0 2 

2.9X10 -12 

0.0 

-320.0 

H0 2 4- O = OH + 0 2 

3.0X10 -11 

0.0 

400.0 

II0 2 + H = OH 4- OH 

4.2 X 10 — 10 

0.0 

1900.0 

H0 2 + H = H 2 + 0 2 

4.2X10 -11 

0.0 

-700.0 

ho 2 + ho 2 = h 2 o 2 + o 2 

2.3X10 -13 

0.0 

1200.0 

h 2 o 2 + 0 = OH + ho 2 

1 .4 X 10 — 12 

0.0 

-4000.0 

N + O 2 = NO + O 

1.5 X 10 1 * 

0.0 

-7200.0 

N + NO = N 2 + O 

3.4X10 -11 

0.0 

0.0 

NO + O + M = N0 2 + M 

4.7X10 -28 

1.5 

0.0 

NO + H0 2 = N0 2 + OH 

3.7 X 10 — 12 

0.0 

480.0 

NO + NO 3 = N0 2 4- NO 2 

1.7X10 - 11 

0.0 

300.0 

N0 2 + O ss NO + 0 2 

6.5X10 -12 

0.0 

240.0 

N0 2 4-0 + M = NO 3 4- M 

8.1 X 10 — 27 

2.0 

0.0 

N0 2 4“ OH + M = HNO 3 + M 

2.2 X10 ” 22 

3.2 

0.0 

N0 2 4- NO 3 + M = N 2 O 54 - M 

9.9 X 10 — 20 

4.3 

0.0 

NO 3 4 - O = N0 2 4* 0 2 

1.0X10 -11 

0.0 

0.0 

O 4- 0 2 4 - M = O 3 4- M 

3.0 X 10 — 28 

2.3 

0.0 

O 4- O 3 = 0 2 4- 0 2 

8.0 X 10 — 12 

0.0 

-4120.0 

H 4 O 3 = OH 4 0 2 

1.4X10 -10 

0.0 

-940.0 

OH 4 - O 3 = H0 2 4 - 0 2 

1 . 6 X 10 -12 

0.0 

-1880.0 

NO 4- O 3 = N0 2 4- 0 2 

2.0 X 10 ~ 12 

0.0 

-2800.0 

N0 2 4- O 3 = NO 3 4- 0 2 

1.2X10 -13 

0.0 

-4900.0 

OH + HNO 3 = H 2 d + NCTJ'”* 

’ 7.2X10 -15 

0.0 

1570.0 

SO 4- 0 2 = S0 2 4- O 

1.4 X 10 — 13 

0.0 

-4550.0 

SO 4- O 3 = S0 2 + 0 2 

4.5X 10 ~ 12 

0.0 

-2340.0 

S0 2 4 OH 4 M = k§5 3 4- 

"7.5X10“ 2 3 

3.3 

0.0 

HSO 3 + 0 2 = H0 2 4- SO 3 

4.0 X 10 “ 13 

0.0 

0.0 

SO 4 N 0 2 = S 0 2 +N"o"“ 

1.4 X I 0“ 11 

0.0 

0.0 

SO 4 - OH = S0 2 4- H 

8.6 X 10 — 1 1 

0.0 

0.0 

S0 2 4 O 4 M = SO 3 4- M 

4.0 X 10 “ 32 

0.0 

- 2000.0 



TABLE 6. 

Plume Centerline N0 X , S0 X and Oxidizer Speciation 
(Mach 2.4; N47; 18.4 km; June 15) 

species mole fractions 



exit plane 500 ft 

1000 ft 

2000 ft 

OH 

1 . 0x10 

5 1.4X10 7 

1.4 X 10 — 8 

1.9X10 -9 



, —6 


ft 

o 3 

0.0 

2 . 1x10 

2.4X10 

2.4X10 

NO 

4.3X10 

~ 5 6 . 2 X 10 -6 

2.0X10 -6 

9.4 X 10 — 7 

N0 2 

4.8X 10 

_6 6.2 X 10 — 7 

2.0 X 10 — 7 

1.0X10 -7 

NO 3 

0.0 

2.1 X10 -11 

4.5X10 -12 

8.1 X10 -13 

N 2 O 5 

0.0 

2.1 X10 -10 

3.2X10 -10 

3.4 X10 — 10 

HNO 3 

0.0 

2.4 X 10 — 7 

1.1 X10 -7 

5.9 X 10 — 8 

S0 2 

6.9X10 

_6 9.5 X 10 — 7 

3.0 X 10 — 7 

1 .4 X 10 — 7 

SO 3 

0.0 

6.8X10 -8 

3.2 X 10 — 8 

1.6 X 10 — 8 

NO 2 /NO 

0.11 

0.10 

0.10 

0.11 

hno 3 /no x 

— 

0.035 

0.05 

0.056 

SO 3 /S 02 

— 

0.07 

0.11 

0.11 


The SO 3 and N 2 Os created in the exhaust plume will react immediately on contact with contrail 
droplets or particles to form condensed phase H 2 S0 4 and HNO 3 , respectively. Exhaust gaseous 
HNO 3 will also condense onto contrail particles. Each of these species, gaseous SO 3 , N 2 Os and 
HNO 3 can also be expected to interact with exhaust soot, as can exhaust N0 2 23 and entrained 
atmospheric O 3, 24 to create oxidized soot surfaces capable of nucleating either exhaust contrails 
or atmospheric acid/water aerosols after plume mix-out. This soot surface conditioning may be 
an important plume chemical effect since only 1% of the soot particles formed by burning jet 
fuel can immediately condense water vapor at typical plume supersaturation ratios without such 
conditioning. 25 We plan to incorporate these heterogeneous processes in future plume chemistry 
models. 


VORTEX WAKE DYNAMICS 


The scales and regimes involved in the dispersion of pollutants from the exhaust of stratospheric 
aircraft have been discussed by several authors mainly in connection with the CIAP studies in the 
70’s. 26,27 The early part of the process (the first hour) was described in some detail by Hoshizaki 
et a i., 28 who called this period “the wake regime". The first few tens of seconds of this period 
are characterized by the dominant role of the aerodynamic flowfield of the aircraft (first through 
the momentum of the jet itself, and later through the effect of the aircraft sinking vortex pair). 
Wake dispersion, after an ill-defined transition, is dominated by environmental effects, such as 
wind shear. 

Based on the analysis of horizontal and vertical photographs of B-52 contrails, Ref. 28 concluded 
that in a first phase, called the “jet regime", the engine plumes grow by ordinary turbulent mixing 
to fill the recirculating vortex pair cell (of dimensions roughly 1.2 x 1.5 times the wing span). The 


attraction to the wing tip vortex cores was noted, but dismissed as a minor effect. After this time 
(~10 s), the effluent remains horizontally confined by the cell, but detrainment from the top as the 
cell descends leads to vertical growth. This “vortex" phase ends at about 100 s with the breakup 
of the vortex pair, following which the wake remnants gradually lose their organized rr 'on and 
ambient effects dominate. The effect of the aircraft was noticed, however, for as much as e hour, 

since different aircraft lead to different dispersion rates on that time scale. This concep model 

(along with the models of Holdeman 29 and Nielsen 30 which ignored wake interactions ^re later 
also applied to the analysis of supersonic YF-12 flight data. The interpretation of t' -jTowth of 
this wake is complicated by the visualization method used, as will be discussed belo 

As will be discussed more fully in later sections, the depth of the pressure am. temperature 
minima induced by the rolled-up wing-tip vortices scale as the square of the Mach number, and 
inversely with the square of the wing aspect ratio. Thus, what was found to be a minor effect 
when examining transonic aircraft may become a very important effect for a low aspect ratio 
supersonic aircraft. The horizontal densitometer scans of B-52 contrails (Hoshizaki et al., Ref. 28) 
already show two peaks at the expected vortex core locations, and many other observations (see 
for instance photographs in Ref. 32) indicate the persistence of these two separate, well-defined 
and non-growing columns of condensation material up to the time of vortex breakdown. 

The extent of plume gas trapping in the low pressure troughs caused by the vortical motion, and 
the degree to which exhaust gas is able to escape and be detrained from the sinking recirculating 
vortex pair cell determine the effectiveness of the vortex confinement and mixing suppression. 
These factors are strongly dependent on the relative spanwise location of the engines and the 
vortex cores, and also on the “attraction" of these cores, as measured by their pressure trough 
depth. The trapped and the detrained portions of the exhaust gas follow very different dilution 
and cooling histories. While the classical Appleman 21 treatment (isobaric mixing of effluent and 
ambient air) should be applicable to the detrained portion, the portion captured by the vortex cores 
may undergo significantly lower temperatures during mixing, due to adiabatic expansion into the 
low pressure of the core. Since, to begin with, stratospheric conditions of interest are typically 
within ±10°C of the Appleman contrail threshold, these deviations may well initiate ice formation 
in the cores under conditions for which contrails would not be predicted ordinarily. In addition, the 
trapped gas dilutes at a much reduced rate, and this may provide ice growth times and/or chemical 
reaction times much longer than otherwise expected. 

Given the stronger confining effect of the vortex cores in a supersonic, as compared to a 
transonic, aircraft, it is not even clear whether substantial detrainment from the descending vortex 
pair cell is to be expected in the supersonic case. If vertical detrainment were either absent or 
strongly reduced, the morphology of the effluent distribution at the end of the vortex wake would 
be quite different from what has been accepted since the studies in the 1970’s. In the YF-12 data 
of Ref. 31, there was little vertical growth in the wake regime, amounting to as little as 1/4 of that 
predicted by the “detrainment” model of Ref. 28. Unfortunately, interpretation of these supersonic 
cruise data is complicated by the fact that the measurement was based on in-flight fuel dumps from 
the aircraft centerline. Only those portions of the dump which mixed with the warm buoyant plume 
would be subject to vortex capture, leaving the rest free to detrain. 

An additional difference between conventional and supersonic aircraft wakes must be mentioned 
here. For most conventional aircraft, the vortex pair becomes unstable and breaks up into irregular 
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vortex rings at about the same time as its vertical motion is arrested by the stable atmospheric 
stratification. Organized vertical motion would in any case be ended by the latter effect at this 
point, and the coincidence has led to an identification of “vortex pair breakup" with “vorticity 
dissipation". In the supersonic case, however, the vortex pair has more vertical momentum, and 
its vertical motion will continue (in the form of vortex rings) well after the straight vortex system 
has broken up. This will probably lead to pollutant deposition at lower altitudes (by a few hundred 
meters) than expected, and also to longer than expected confinement times. The effect may also be 
significant in terms of wake hazards to other aircraft. 

The following section, plus Appendices A and B, provide a more detailed discussion of these 
effects, based on a series of simplified analytical models which are intended to approximately 
quantify their relative importance and scales. 

Analysis 

Aircraft and Vortex Parameters Table 7 summarizes the parameters assumed for our 
vortex wake analysis. The airplane is a Mach number 2.4 design beginning its cruise at 17.4 km 
altitude. For estimates of wake properties, an elliptic loading distribution has been adopted, and 
the vortex roll-up distance was estimated based on Ref. 33. 


TABLE 7. 

Assumed Aircraft Data for Wake Dynamics 
2 

Wing area: A \y — 660 m Span: b = 39 in 

Weight (cruise): W = 2.76 X 10^ N 

L 2 

Flight altitude: Z — 17.4 km Aspect Ratio: "Aw ~ ^ 

Mach No.: M = 2.4 Speed: V = 708 m/s 

Atmosphere: p — 0.135 kg/m^ f T — 217 K , p — 83.8 mb 

Lift coefficient: C L = 0.124 

Exhaust temperature: T e = 445 K 

Fuel/air ratio; 0.0155 

Airflow: m a = 97 Kg/s/engine 

Exhaust speed (relative): 1430 m/s 

Wake time scale: T = = ^~W ~ = 2.04 s 

Vorticity roll-up time: tr.u. — 1.5 T = 3.07 s 

Vorticity roll-up distance: X ru . = 2170 m (56 b) 

Rolled-up vortex spacing: 2 y = 'jt = 30.6 m 

Centerline vortex circulation: To = j^=943 m^/s 

r 

Self-induced descent velocity: Wo = 2n(2y) = “*- 90 m / s 

Vortex break-up time: tv.b. — 15t = 30.7 s 

Vortex break-up distance: %v.b. ~ >.b. = 21,700 m (556 b) 


The structure of the vortices after roll-up has been calculated based on the Betz-Donaldson 
model, 33,34 according to which the circulation originally bound to the wing outboard of a spanwise 
location t/i will end up rolled into the inner part of the trailing vortex core, within a radius r from its 
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center — which equals the distance from y\ to the centroid of the circulation distribution outboard 
of yi . For the assumed elliptic loading, this model predicts a tangential velocity distribution near 
each core given by 



where v is the flight speed, Cl the lift coefficient, AR the wing aspect ratio, b the wing span and r 
the distance to the vortex core. This has been confirmed in tests, 33 and it appears to hold to radii 
much smaller than the classical core size estimate of Prandtl, 34 with the result that v reaches higher 
values near the core than expected on the basis of a model with all the vorticity concentrated in the 
Prandtl core. The validity of the Betz-Donaldson model has been debated, 35 and no rigorous basis 
for it seems to exist. We use it here simply as a convenient semi-empirical device. 

Using Eq. (1), several quantities of interest can be calculated. The centripetal acceleration vj/r 
is, using the parameters in Table 7, 


y centripetal — ^2 0^ MKS Units) . (2) 

With reference to Fig. 1, this gives 26 m/s 2 at the inboard plume location, and 73 m/s 2 at 
the outboard plume location. This shows that centripetal buoyancy effects will be strong and 
will compete with (perhaps dominate) those of upwards buoyancy and turbulent detrainment of 
effluents from the cell’s edge. Notice (Eq. (1)) that the peculiar features of the HSCT (small AR, 
high u) greatly emphasize this effect. For a transonic airplane of the same wing area but with a lift 
coefficient of 0.5 and an AR of 10, g centripetal would be less by a factor of 17.2 at the same value 
of b/r. This makes plume/wake interaction a secondary issue for such aircraft, and may account 
for its dismissal by Hoshizaki et a i. 28 and, more generally, for the lack of attention it has elicited 
so far. 

The pressure gradient set up by the vortex follows easily from the hydrostatic radial force 
balance dp/dr = —pvj/r. From Eq. (1) and ignoring the relatively small effects of compressibility 
the local depression, A p = Poo - P. is seen to vary as 1/r. Assuming adiabatic expansion, the 
corresponding temperature depression is found to be 

AT 6(t - 1) M*Cl b 

7„ ir< AR? r ' ' 1 


6.9m 1 1 .5 m 



Fig. 1 . Engine location with respect to rolled-up vortex core. The core is at (irf 4)(b/2) 
from centerline. 
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For the parameters of Table 7, we find A T(K) = -3.5 jr{m). Once again, for typical transonic 
aircraft the corresponding AT is less by a factor of 83, making it almost unobservable. 


The vortex system is known to break up due to growth of the Crow instability (vortex bursting is 
unlikely at the moderate C L of cruise). The estimate for breakup time given in Table 7 {t B = 15r) 
is consistent with data of Ref. 36 for a quiescent atmosphere, and is 50% longer than the estimate 
given by Hoshizaki et al. 28 for tropospheric conditions. 


Turbulent Plume Growth with Buoyancy Effects The initial spread of the engine 
plumes must be controlled by turbulence generated by the plume/ambient velocity shear. However, 
centripetal buoyancy effects must be felt strongly, as noted, and will eventually dominate the plume 
dynamics. The overall dynamics is fairly complex, and we will attack it initially by considering 
simpler model problems. The first question is the estimation of the time after which buoyancy 
dominates over jet growth. To this end, we temporarily ignore the shearing of the plume by the 
vortex flowfield, and consider a buoyant round engine jet in the co-flowing airstream. The analysis 
is relegated to Appendix A, and is based on the use of a turbulent diffusivity model which combines 
those of Prandtl 37 for a coflowing jet and of Morton, Taylor and Turner 38 for a rising linear plume. 
The resulting jet diameter is found to evolve with distance x as. 


where 
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[l+0G(Ol‘ /3 


£ = x/e 
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1 T]ov 9eff 


d = 0.061 





(4a) 

(46) 

(4c) 


and 

G(0 = tVuT 2 +ln(Z + \/l + £ 2 ) • (4<0 


Here, rj ov is the overall engine efficiency (fuel heat value to thrust), F is the engine thrust and 
Do is the engine nozzle diameter. The acceleration, g e /f, could be either the actual gravitational 
acceleration, or the centripetal acceleration calculated above. 

By examination of Eqs. (4a)-(4d) we can see that when ( « 1, 9 ~ x l D (jet-dominated 
regime), and when £ > 1, D = x 2 / 3 (buoyant plume regime). The two are nearly equivalent at 
£ = 2. For the parameters of Table 7, and using the centripetal g for either engine, we find that 
£ = 2 is reached when D/D 0 =2.17 (outboard engine) or D fD 0 = 2.89 (inboard engine). Thus, 
by the time centripetal buoyancy becomes the dominant plume growth mechanism, the dilution 
ratio (proportional to ( D/D 0 ) 2 ) is about 4.7 to 8.3, depending on which engine is considered. 
For comparison, using “natural" gravity (g =9.8 m/s 2 ), the crossover occurs when D/D 0 = 4.60 
(dilution ratio = 21), by which time the plumes would fill a substantial fraction of the vortex cell. 

The analysis above has not accounted for the effects of compressibility on mrbulent jet mixing. 
Ref. 6 shows substantial mixing suppression when the “convective Mach number", M* exceeds 
about 0.5. M* is measured with respect to an intermediate frame in which both pressure and total 
pressure are equal for both streams. In the case of equal specific heat ratios, this reduces to 


M* = 


v? ~ vi 

02 + Ol 


(5) 
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where v and a are the respective flow and sonic velocities of the two streams. In our case (Table 7) 
M* = 1.01, and the initial mixing rate should be reduced by a factor f c between 0.2 and O.4. 6 The 
factor f c will then approach unity as mixing proceeds and the relative Mach number decreases. 
The effect on the above analysis would be that the average f c would multiply both 6 and i. 

A Model for Plume- Vortex Interactions The fluid mechanics of the turbulent buoy- 
ant jet in the flow field of a concentrated vortex is too complex to be adequately analyzed here, 
but we will at least account for the salient features in an attempt to estimate the capture time of the 
engine exhaust into the vortex cores. The following assumptions are made: 

a. The plume is convected and sheared by the (un-modified) Betz-Donaldson vortex flow. The 
shearing spreads the plume around a circle centered at the vortex core, at a rate proportional 
to the difference v 9 (r - D/2) - vg(r + D/2) where r is the distance to the core and D is 
the transverse plume size (See Fig. 2). 

b. The plume cross-section grows by entrainment of ambient air, in proportion to its inward 
radial velocity times its projected area (Morton-Taylor-Tumer model 38 ). No account is 
taken of axial flow, except that the model is used to estimate a plume diameter as it begins 
to be affected by the vortex. 

c. The radial motion is obtained from a momentum balance, using the centripetal buoyancy 
ofEq. 2. 



Fig. 2. Geometry of plume trapping in a plane transverse to the flight path. 

Letting A be the plume cross-section at a given time, and A p the mean density defect in the 


plume, mass conservation imposes 

o 

II 

<1 

(6) 

The entrainment rate assumption gives 

dA 

( 7 ) 


— = a As iv . 
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where As is the tangential plume spread (Fig. 2), v r is its radial velocity (inwards), and a = 0.94 
is chosen to correspond for the circular cylinder case to the value used in Ref. 28. 

Equating the buoyancy force to the rate of momentum increase due to mass accretion gives 


otp a v 2 T As = — AAp . (8) 

r 

Finally, the shearing assumption is expressed as 



where D is assumed constant. 

Combining Eqs. (6)-(9) with the expression Eq. ( 1 ) for vg leads to a set of differential equations, 
for which closed-form solutions can be obtained: 


and 
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( 11 ) 

( 12 ) 

(13) 

(14) 


For application, the initial distance r 0 is 6.9 m for the outboard and 11.5m for the inboard engine. 
The initial diameter D 0 and dilution (or density depression {Ap/ p)o) can only be approximately 
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selected to account for the turbulent growth rate in the initial jet stage. A rough estimate, based on 
the arguments of turbulent plume growth (ignoring compressibility effects) gives (for an engine jet 
diameter of 2 m) initial diameters of 4 m and 5 m for the outboard and inboard engine respectively. 
Correspondingly, the density depressions (Ap/p)o are 0.122 and 0.0816, since the value at the 
engine face is 0.51. 

Figs. 3 and 4 show the subsequent evolution of plume-core distance r, angular spread Ad of 
the sheared plume and overall dilution ratio A/Aj et . Since the radial pressure gradient intensifies 
towards the vortex center, the radial velocity of the plume is seen to increase with time, despite 
the progressive dilution. Thus, although it is not clear what distance r to pick as the end of the 
“capture" process, one can make a reasonable estimate of the “capture time". For instance, Eq. (10) 
indicates a finite time for r — ► 0, and this can also be seen from the downwards turning of the r(t) 
curves in Figs. 3 and 4. Perhaps a more reasonable capture time can be chosen as that at which r 
equals the original plume diameter D 0 (after allowing for the initial turbulent spread) or, what is 
nearly equivalent, that at which the plume has been fully spread into an annulus (Ad = 360°). This 
gives for the outboard engine about 90 wing spans, by which time the overall dilution (including 
initial turbulent growth) is about A/Aj et = 28. The process is slower for the inboard engine, and 
it takes 350 wingspans to capture that plume, with an overall dilution ratio of 56. 



Fig. 3. 


Inner engine jet mixing and displacement. Initial jet diameter, 
distance from vortex core, r q = 11.5 m, and density difference, 
0.0816. 


Do = 5 m, 
(Ap/p) o = 
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Fig. 4. Outer engine jet mixing and displacement. Initial jet diameter, D 0 = 4 m, 
distance from vortex core, r 0 = 6.9 m, and density difference, (A p/p)o = 0.122. 

The basic wake time scale (Table 7) is r = (2 AR/C L )(b/V), in terms of which roll-up 
takes approximately 1.5r and break-up takes 15r (or, for our values of AR and Cl, 56 and 560 
wingspans, respectively). We see therefore that a) the time for the plume capture process scales 
also as the basic time r (see Eq. (10)), and b) for both engines, the capture time is intermediate 
between those for roll-up and break-up, as required for validity of the modeling approach. 

It is interesting to observe, in connection with point a) above, that, for a given geometry and 
initial density defect, plume capture occurs at a fixed fraction of wake breakup time, independent of 
Mach number. Thus, of the two factors which are clearly different between transonic and supersonic 
aircraft (AR and M), only one, .47?, remains to shorten the capture time in relation to wake lifetime. 
The reason why flight speed does not influence the time ratio is that all processes, including buoyant 
radial velocity and azimuthal vortex velocity, scale up together as v (see Eqs. (1) and (8)). This 
fact can be exploited for setting up experimental simulations of these phenomena in water tanks. 
The Mach number is still important however, in dictating deeper temperature and pressure minima 
in the supersonic case, and hence strengthening centripetal buoyancy when compared to natural 
buoyancy. Whether or not these stronger temperature depressions lead to enhanced condensation 
depends on the degree of dilution ultimately achieved by the captured plume. 
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The Last Stages of the Vortex Regime 

As noted in Table 7, the wake is expected to survive as an organized cylindrical vortex cell 
for about 22 s. According to our analysis in above, the vortex cores at breakup time may contain 
a large fraction of the engine exhaust, at a dilution smaller than would have been expected in the 
absence of the vortices, and at a few degrees lower temperature. In attempting to extrapolate from 
this point, several questions arise: 

a. For how much longer will the confining effects of vorticity prevent or delay plume disper- 
sion? 

b. Under what conditions will the higher plume vapor concentration and lower temperature 
lead to contrail formation? 

c. If condensation does occur, will the confinement time be long enough for significant 
heterogeneous chemical interactions to occur or for particles to grow to precipitation sizes? 


Vortex Dynamics After Break-up Unfortunately, the fluid mechanics of the broken- 
up vortex system is still poorly understood, and so, only very rough estimates can be made on the 
important question a). Experimental and analytical stability results 36 indicate that reconnection 
due to growth of the Crow instability tends to produce irregular elongated ring vortices with aspect 
ratios of the order of 4:1. Dhanak and DeBemardinis 39 followed numerically the evolution of an 
elliptical ring vortex of this type and found that, after an oscillation in which the short and long 
axes interchange positions, the cores touch in the middle and reconnection can be expected, leading 
to two roughly 2: 1 rings. These are then structurally stable, and oscillate as described. The time 
for subdivision of the original ring is (for our parameters) about 5 s. 

If we can ignore temporarily the disturbing effects of wind shear, the smaller rings may continue 
to descend for some time, as they entrain new air by turbulent diffusion, and consequently grow 
in size and slow down. Glezer and Coles 40 did careful measurements of the evolution of circular 
turbulent rings in water, and showed that similarity exists in this motion provided a “virtual origin" 
is identified. In Appendix B we show that their results can be explained by means of a simple 
entrainment hypothesis, for which the entrainment parameter is derived from the data. As explained 
in Appendix B, the linear dimensions of the ring increase as (1 + t/<i) 1//4 , where t\ = nao /2/3 vq 
(a 0 = initial core radius, v 0 = initial ring velocity, and /? = 0.01). For an estimation, assume in 
our case an initial ring diameter 1.5 times the linear vortex spacing (i.e., 48 m), and an initial core 
diameter 0.1 of the ring diameter (a 0 = 2.4 m). Since v 0 must be close to the linear vortex pair 
velocity of 4.9 m/s, we calculate t\ =76 s. 

The effects of atmospheric stratification must be also considered at this point. As the air mass 
enclosed within the vortex cell sinks into a stably stratified atmosphere, it will undergo adiabatic 
warming due to the higher pressures encountered, and will develop buoyancy. The effect of this 
buoyancy will be to eventually lead to an oscillatory vertical motion at the Brunt-VSisala frequency, 


N = 



(15) 
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( 16 ) 


where ( -dT a /dz ) is the atmospheric lapse rate, and 

_dT _ 7-1 T_ 
dz 7 H p 

is the warming rate of the descending air (H p is the pressure scale height of the atmosphere). 
However, a second effect of the descent-induced buoyancy is of interest to us. This is the production 
of vorticity opposite in sense to that of the wing vortices. The mechanism is clear from the 2 - 
dimensional vorticity equation: 

M£ + S-V(-) =-i V(-) xVp . (17) 

at p p p 

The pressure gradient Vp points vertically down, while -Vl/p is concentrated on the cell’s 
edges, and points outwards from the cell. Thus, baroclinic vorticity is generating along the cell’s 
edges, in the sense contrary to that in the vortex cores. At the time buoyancy stops the cell’s descent, 
the total counter-vorticity created is also sufficient to cancel that of the wing vortex system, although 
its distribution is different. We can therefore expect the vortex system, whether in the form of the 
undisturbed cylindrical vortex cell, or of its successor ring vortices, to vanish in the vicinity of the 
lowest point of the Brunt- VSisala cycle. This occurs when Nt = n/2, with the vertical descent 
distance being then A z = -w 0 /N (w 0 = initial descent velocity). This behavior can be clearly 
seen in the water-tank data of Ref. 36, and was also assumed by, for example, Greene 41 in his 
model of wake decay. 

For our parameters, with dT a jdz = 0, H p = 6380 m, we calculate N = 0.0210 radians/s, and 
the stratification-induced vortex destruction can be expected to happen at t = ir/2N = 75 s, after 
a total descent of — Az = 232 m. Notice that this is 2.4 times longer than the time required for 
Crow instability to occur, so that a substantial period (44 s) can be anticipated during which the ring 
vortices formed after breakup can retain their individuality and continue to trap the engine effluent. 
The amount of turbulent entrainment into the rings in this time can now be estimated from our 
Appendix B results. With 1 1 = 76 s, the ring linear dimensions will increase by (1 + 44/76) 1 / 4 = 
1.12, with additional dilution by a factor (1.1 2 ) 3 = 1.41. 

Condensation Considerations A slight modification of the classical Appleman argu- 
ment 21 concerning contrail formation will help illustrate the potential vortex effects. Assuming 
that each kg of burnt fuel produces 1.29 kg of water vapor and 43 MJ of heat, of which the fraction 
1 — rj overa u = 0.5 appears as sensible heat in the plume, and that the plume dilutes continuously 
by isobarically mixing with air at temperature T a and moisture w a (g/kg), its average temperature 
and moisture will evolve together according to 

T = T a + 17.7(w — w a ) . (18) 

In Fig. 5, we have represented in a T-w plot the water and ice vapor saturation lines (for an air 
density of 0.135 kg/m 3 ). The plume material will, according to Eq. (18) evolve down an inclined 
straight line of slope 1/17.7 towards the existing atmospheric conditions (T a , w a ). If this line 
intersects the water saturation line, droplets will condense on existing condensation nuclei, and 
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immediately freeze to form ice particles. These then persist even as the plume dilutes below the 
water saturation line, and will re-evaporate only after the ice saturation line is re-crossed. 

In the presence of the vortex depression, however, both the ambient air and the ime gas are 
expanding and adiabatically cooling as they mix. Ignoring differences in the sped' heat ratios, 
Eq. (18) must then be modified to 


T = 


T a + 17.7(U> - W a ) 



(19) 


where 7 = 1.4. Tracing the plume trajectory in the T-w plane now requires knowledge of the 
relationship between local pressure at the plume location and degree of dilution at that point. 



Fig. 5. Temperature — Water vapor content (T-w) plane with water and ice s ration 
curves overlaying exhaust dilution trajectories. 

As an example, consider a dry stratosphere at T a = 215 K which, as Fig. 5 she is at the 
threshold for contrail formation for this air density (plume expansion along abc). Usir, le p me- 
vortex interaction model we can calculate at each core-plume distance r the correspom * dilution 
ratio A/A 0 and temperature depression (Eq. (3)). Starting from an engine-exit val w = 20 
g/kg, we can then calculate the new plume trajectory a'b'c (Fig. 5), which shows defii e contrail 
formation. The last portion of the a'b'c curve is shown dotted, as it corresponds to ihe poorly 
understood process of vortex dissipation and final dilution. 
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The above results indicate an increase by about 2°C of the minimum atmospheric temperature 
required to form a contrail. The possible global significance of this is best appreciated by reference 
to Fig. 6, which is the result of straightforward application of Appleman’s criterion. The region 
labelled ALWAYS corresponds to temperatures lower than that at points like c in Fig. 5, where 
contrails form even in dry air. The region labelled NEVER corresponds to temperatures greater 
than T b in Fig. 5, where contrails would not form even in a saturated atmosphere. By inspection 
of Fig. 6, the cruising altitude range (17-20 km) for a Mach 2.4 aircraft contains a significant 
proportion of conditions lying no more than a few degrees from the ALWAYS limit. Thus there 
could be a noticeable impact on the frequency of contrail formation associated with wake-induced 
cooling. It may be, however, that a more significant effect than this contrail probability increase is 
the reduced dilution and extra confinement time provided by the wake vortices. 



March 22 



Fig. 6. Predictions using equilibrium, fully-mixed condensation algorithm. 

The size to which the contrail ice particles can grow is mainly determined by the concentration 
of active condensation nuclei (cn) in the plume; kinetic limitations are less important, except 
perhaps for diameters over 10 p. Unfortunately, not much is known about cn production by 
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engines, and even less about their activation mechanisms. Knollenberg 42 sampled a Sabreliner 
contrail in saturated ice conditions and estimated 10 10 /m -3 active cns at the engine exhaust plane. 
Rosen and McGreegor 43 measured photoelectrically the total carbon particle concentration some 
600 m behind an F-104 plane. Their data seem to indicate a total particle concentration of about 
2.5 x 10 lo /m -3 extrapolated to the engine exhaust. On the other hand, Hoffmann and R >sen 44 
encountered an 18 hour old SR-71 contrail which had sheared to a horizontal width of 2C 0 km 
(but with a vertical thickness of only 200 m), and estimated emission of 3 x 10 13 to 12 10 14 

m -3 particles greater than 0.01 pm (not necessarily all active). The fraction of all particles emitted 
which become active cns is also uncertain. Hailed et a l. 25 measured a conversion fraction of 1 to 
3% for JP-4 fuel, but the 3% level was only achieved after aging for 20 hours. 

From the example of Fig. 6, it is possible to have about 50% of the original vapor in the form 
of ice particles when the ambient temperature is near the contrail threshold. Assuming also 1.29 
kg of water/kg of fuel and p = 50 mb, T ex it = 450 K, and letting n (m -3 ) be the active nucleus 
concentration at the engine exhaust, we estimate an eventual mean particle size 

m ) = (^)‘ /3 ( 20 ) 

which, for the range of estimates of n given above implies R values between 2 and 0.06 pm. A 2 pm 
radius particle of pure water, falling through dry air at 210 K, would evaporate fully in about 10 s 
and settle no more than about 1 cm. This would indicate the lack of any contrail settling. However 
two modifying factors may occur in the case of an HSCT contrail, a) Lean-burner combustors, 
being developed for emissions reduction, will produce much less soot than conventional burners — 
perhaps by two or more orders of magnitude. Active ccn levels may be significantly lower still 
Hallett et ai. 25 . This may lead to ccn concentrations approaching background levels and hence 
(Eq. 20) to much larger particles, b) The extended confinement of exhaust products may lead to 
the formation of significant surface layers of nitric acid trihydrate (NAT) and other hydrated acid 
species on the ice particles (see Eqs. 2 1 -28). This would lower their vapor pressure and essentially 
prevent re-evaporation, with the result of a much larger settling distance for a given particle size. 

The acid gas processing of the condensed exhaust water begins with the plume/wake chemical 
effects described in the previous section. The chemistry modeled there includes the formation of 
acid gases and their precursors (HNO3/N2O5 and H 2 S 0 4 /S 0 3 ). These gases are formed through 
the scavenging of exhaust OH radicals and reaction with entrained atmospheric O3. They are 
then available for heterogeneous chemistry processes that can change exhaust speciation through 
processes like: 


N 2 0 5 + H?0 (c) — ♦ HN0 3 (c) , 

( 21 ) 

S0 3 + H 2 0 (c) — * H 2 SO 4 (c) , 

( 22 ) 

HN0 3 + H 2 0 (c) -> HN0 3 (c) , 

(23) 

H2SO4 + H 2 0 (c) -» H2SO4 (c) , 

(24) 

HC1 + H 2 0 (c) -+ HC1 (c) ; 

(25) 


temporarily removing plume and atmospheric acid gases from the gas phase or setting up the 
liberation of reservoir atmospheric chlorine: 

CIONO2 + H 2 0 (c) HOC1 (c,g) + HNO3 (c) , (26) 
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( 27 ) 

( 28 ) 


C10N0 2 + HC1 (c) -> Cl 2 (g) + HN0 3 (c) , 

HC1 + HOC1 (c) Cl 2 (g) + H 2 0 (c) . 

The activation of plume/exhaust soot to ccn will occur on particles through their reaction with 
acid and oxidizing gases. In future work, we will model both this soot activation process and 
the kinetic nucleation and growth of contrail particles, a process which depends directly on the 
degree of plume/wake water vapor supersaturation, determined by the plume/wake temperature and 
entrainment/detrainment profiles and the number of active condensation nuclei. Particle formation 
and growth through the wake dispersal regime depend on more than the simple equilibrium H 2 0 
condensation considered in the plume mixing and chemistry section. Clearly, in this heterogeneous 
chemical environment, the relevant condensation phenomena determine the degree to which gravi- 
tational settling can drop the particle and its associated acid gases and exhaust particulates to lower 
altitudes less conducive to ozone depletion. 

SUMMARY AND CONCLUSIONS 

The mixing, chemistry and condensation in the exhaust flowfield behind an HSCT were modeled 
as a first step in following the engine emissions to their eventual deposition in the atmosphere. 
Finite rate chemistry was simulated internal to an engine for conditions that are repre sentative of 
engine cycles that would be used in propelling an HSCT. This simulation was used to calculate 
exhaust emissions at the engine exit plane where super-equilibrium concentrations of CO, NO, 
N 0 2 , HNO3 and OH were estimated to be present. , 

The exhaust emissions were followed from the engine exit plane using the standard plume 
flowfield code, SPF-2. The plume code was initialized using estimated emission indices or calcu- 
lated exit plane concentrations. The chemical kinetics in the plume code were upgraded for these 
calculations to include the chemistry for atmospherically active species. N 0 X in the plume was 
partially convened to HNO3 through reaction with the OH radical to leave over 5% of the NO x as 
HNO3 600 m downstream of the exhaust exit. S 0 2 was also oxidized (to SO3) so that in excess of 
10% of the SO x is available for direct conversion to H 2 SC>4 upon contact with condensed H 2 0 . 

Equilibrium water condensation was also added to the upgraded plume code to begin to include 
some of the condensation phenomena needed to incorporate heterogeneous chemical reactivity in 
the plume. Contrail formation was predicted for several altitude-latitude-season combinations. 
These predictions based on equilibrium H 2 0 condensation were consistent with the Appleman 
algorithm for contrail occurrence. 

The SO3 , N 2 C>5 and HNO3 produced in the plume react and/or condense on the aqueous exhaust 
aerosols to form condensed acid solutions. These species as well as N0 2 and entrained atmospheric 
O3 can also oxidize carbonaceous soot particles emitted with the exhaust. The oxidative activation 
of additional nucleation sites affects the size and quantity of condensed particles available for 
heterogeneous chemistry. The aerosols formed and activated in the plume are subjected to the 
vortical flowfield in the wake of the airplane, which transports and processes the particles and gases 
until the wake disperses. - 1 (■■■ • • .• •• 

The vortex wake due to a Mach 2.4 airplane was analyzed to estimate its effects on the transport 
and mixing of the exhaust emissions. For parameters representative of a low aspect ratio supersonic 
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aircraft, the circulation due to the rolled-up vorticity convects and shears the exhaust plume as it 
wraps the exhaust gases around the core. Since the emissions are warmer than the ambient, the 
pressure gradient generated by the vorticity attracts the less dense gas toward the core. 

The vorticity has rolled up into a concentrated trailing vortex pair within about 56 wingspans 
behind the airplane. This vortex wake then propagates in the stratified atmosphere until the mutual 
interaction between the two vortices of the pair cause reconnection and the formation of elongated 
elliptical vortex rings at about 560 wingspans. The “capture" of the exhaust plumes by the vortices 
is dependent on engine placement but, for the case considered, the outboard exhaust is captured 
at about 90 wingspans while the inboard exhaust takes 350 wingspans. These values depend on 
the airplane configuration (aspect ratio) but not on its speed (Mach number). However, the speed 
does affect the intensity of the vortices, their pressure field, and thus the attraction and possibly the 
confinement of the exhaust. 

The same chemistry that occurs in the plume continues as the exhaust is transported and mixed 
by the vorticity in the wake until wake dispersal, with subsequent atmospheric transport. The 
kinetics used in calculating the evolution of the plume will be incorporated into the vortex wake 
model to follow the chemistry from the capture of the exhaust through the dispersal of the wake. The 
simplified analysis performed in the present study provides insight and estimates of the dominant 
processes that are occurring throughout the wake and further can be used to generate first order 
estimates of the flow properties. However it is apparent that the several individual vorticity-induced 
phenomena are, in fact, strongly coupled and a more detailed, presumably numerical, analysis will 
be necessary for quantitative estimation of the evolution of the emitted exhaust gases and their 
transport. 
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APPENDIX A 


Buoyant Turbulent Jet in a Co-flow 


Consider the jet depicted in Fig. Al, issuing with some velocity and temperature excess into a 
parallel stream. A top-hat model of the distributions will be adopted for simplicity. The fluxes of 
axial momentum and of enthalpy will be conserved, while the vertical momentum flux increases 
with distance due to the buoyancy. The convection velocity will be approximated by the external 
velocity u e . We then have 


p e U e [u - U e ) 


7 VD 2 

4 


= F = constant 


(A.l) 
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These equations must be supplemented by one which describes the turbulent diffusion and its effect 
on jet growth. In general, 

u e ^=4 D, (.4.4) 

ax 
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Fig . A 1 . Buoyant jet geometry. 

Where D t is the “turbulent diffusivity". For a non-buoyant jet, a reasonable approximation 
is 37 D t = kD{u - u e ), with k ~ 0.016. For a non-flowing cylindrical plume the Morton-Taylor- 
Tumer 38 model gives D t = aDw (a ~ 0.3). These two limiting forms can be interpolated by 
assuming 

D t = I\D y/a 2 w 2 + (it — u e ) 2 , (A. 5) 


where a = a/K ~ 19. 

For the case of an engine exhaust in flight (u e = v), F represents the engine’s thrust, and Q its 
thermal energy output. Hence, if r] ov is the overall propulsive efficiency. 


F U t T) ov 

Q 1 - Vov 


(A.6) 


To integrate these equations, u — u e from Eq. (A.l) and w from Eq. (A. 3) (after integrating) 
are substituted into (A.5) to give D t as a function of D and x. This is then substituted in (A.4), 
and integration gives the results quoted in Eqs. (4a)-(4d). 


APPENDIX B 

A Model for Turbulent Vortex Ring Evolution 

The similarity results obtained experimentally by Glezer and Coles 40 for the growth and slowing 
down of a turbulent vortex ring can be understood and extrapolated by a simple model based on an 
entrainment hypothesis. We assume: 


a. The vortex core cross section increases at a rate proportional to its own radius and to the 
ring velocity. With reference to Fig. Bl., 

d(T:a 2 ) . . 

dt =/3aV • (£U) 

where /? is an empirical constant. 

b. Geometrical similarity is preserved, i.e., R and a grow in the same proportion: 


a R 

oq Ro 


(B. 2 ) 


c. The vertical momentum is preserved. In combination with (B.2), this gives 

a?v — OqUo ■ (5.3) 



Fig . B 1 . Geometry of descending vortex ring. 

Eqs. (B.l), (B.2) and (B.3) can be combined to obtain explicit laws for the time histories of the 
various quantities, including the distance x travelled by the ring. Ref. 37 introduced the distance 
x n and time to which act as the origin for the measured self-similarity. These quantities are both 
negative, i.e., the “virtual origin" occurs upstream of the point of creation of the vortex ring. Our 
analytical results, of the same form as the data trends, are 


where 


a R , t ,i/4 

ao Ro to 
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Vo to 
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In Glezer’s experiments, 40 v 0 = 128 cm/s for the piston producing the rings. For the ring itself, 
we take v 0 = 64 m/s. Also, R 0 = 0.953 cm, and a 0 ~ 0.2/2 0 , as estimated from the data. The 
measured virtual origin parameters were 

t 0 = -0.44 s , x 0 = -145 cm . 

Using these results, (B.5) gives 0 = 0.00825, while (B.6) gives 0 = 0.0106. Thus the model 
appears to give a reasonable description of the data using 

0 ~ 0.01 . (5.9) 
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Chapter 3 

Plume Chemistry and Dispersion Modeling 
to Evaluate the 

Atmospheric Effects of Stratospheric Aircraft * 

INTRODUCTION 

Our modeling of the dispersion of the wake has concentrated on: a) developing, testing and 
utilizing a Lagrangian air parcel model to simulate the behavior of chemical constituents in the wake 
regime of a plume; b) modifying and utilizing an existing radiation code at AER to calculate the net 
heating/cooling of aircraft plumes in the same regime. The wake regime starts after breakdown of 
the plume vortices; expansion of the plume is still controlled by internal energy dissipation in the 
wake regime. After about 1,000 - 3,000 seconds, the plume expansion and motion are controlled 
primarily by the dynamics of the background atmosphere. 

Lagrangian Box Model 

The Lagrangian box model solves an equation of the form (Gelinas and Walton, 1974): 

= * -Li- ( dCn j' W ) h(«) - >*?(<)] + R<(‘) , (i) 

where Pi and L t denote the chemical production and loss for species i, rii(t ) and n“(f) represent 
the concentrations of species i within the plume and in the background atmosphere, respectively, 
and R, denotes other production/loss processes, such as, for example, denitrification. The factor 
din V(t)/dt represents the entrainment parameter, where the volume V(t) is defined as that 
physical volume which contains all the initial plume molecules. 

Previous calculations of the photochemistry in the wake regime (Gelinas and Walton, 1974; 
Hoshizaki et ad., 1975) concluded that very little photochemical activity took place during the wake 
regime, and the time development of the mixing ratios of trace species was controlled primarily by 
plume expansion and entrainment of air from the background atmosphere. Chemical schemes have 
however changed since that period of time, particularly in the kinetics of species such as CINO3 
and N 2 O 5 . Furthermore, the recent studies of the photochemistry of the polar stratosphere have 
indicated the importance of heterogeneous reactions occurring in the surface of Polar Stratospheric 
Clouds (PSCs). The following reactions have been suggested (Wofsy, 1978; Solomon et ad., 1986; 
McElroy et ad.. 1986; Tolbert et ai., 1987; 1988a, b; Molina et ad., 1987; Leu, 1988a, b; Rodriguez 
et a I, 1988): 

CINO 3 + H 2 0 -» HOC 1 + HNO 3 (solid, gas) ( 2 ) 

CINO 3 + HC 1 -> Cl 2 + HN0 3 (solid, gas) (3) 

N 2 O 5 + H 2 0 -* 2 HNO 3 (solid, gas) (4) 

N 2 O 5 + HC1 -> CINO 2 + HN0 3 (solid, gas) . (5) 

Reactions (2) - (5) have been observed on ice or nitric acid tri-hydrate particles at very fast 
rates. Fast rates for reaction (4) have been observed on sulfuric acid solutions, with the magnitude 
fairly independent of water content (Mozurkewich and Calvert, 1988; Van Doren et ai., 1991). 


* Work performed by Atmospheric and Environmental Research, Inc. under subcontract to ARI. 
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The new kinetic data could change the conclusions of the assessment undertaken in the mid- 
70’s. In particular, heterogeneous reactions occurring on contrails could alter the partitioning of 
nitrogen species, and favor production of nitric acid. Enhanced concentrations of nitric acid could 
in turn favor production of nitric acid tri-hydrate contrails, particularly since these tri-hydrates form 
at about 5K above the freezing point of water. 

We have developed a Lagrangian box model which includes entrainment of background air. 
This model includes complete gas-phase chemistry of HO x , NO x , Cl x , Br x , ozone, methane and 
ethane. We also include in the model reactions (2) - (5). As discussed below, other heterogeneous 
processes have also been included for purposes of exploration of their impact. The model has been 
used to revisit the issues addressed by the CIAP panel during the mid-70’s. In particular, we have 
addressed and will continue to address the following issues: 

a. What is the partitioning of nitrogen species during the wake regime, and in particular, how 
much NO x is converted to HNO3? 

b. What is the impact of heterogeneous chemistry on the conversion of NO x to HNO3? 

c. What is the impact of entrainment on the chemical development of species within the wake? 

d. What other heterogeneous processes could affect the chemical behavior of species in the 
wake? 

e. What is the ozone loss in the wake? 

f. What is the time development of other radiatively active species, such as water? 

We have carried out exploratory calculations to address some of the above issues. These 
calculations and results are presented in the following. In particular, we look at two regimes: a) 
calculations are carried out for the first hour of the wake regime, assuming very fast entrainment, and 
testing the sensitivity to different assumptions about heterogeneous chemistry; b) the sensitivity 
of the calculations to the entrainment parameterization and scale-dependent dispersion is also 
examined; in particular, we extend some of the calculations to a time period of two days to 
illustrate the time development of important trace species for different dispersion rates. Some of 
the results described below were presented at the HSRP workshop in Virginia Beach, January, 
1991; other results will be presented at the HSRP workshop in Virginia Beach, May, 1992. 

Radiative Calculations 

The final altitude of deposition of aircraft effluents depends on a whole suite of dynamical and 
radiative processes acting at both wake and synoptic scales. The ascent/descent of the aircraft 
plume and its dispersion is determined in part by processes such as buoyant rise, gravitational 
collapse, wind shear, and turbulent dispersion (Hoshizaki et a 1., 1975; Miake-Lye et a 1., 1991). 
The interaction among these processes is still poorly quantified, in part due to the scarcity of data 
on aircraft wakes at stratospheric altitudes. 

Another potentially important process was suggested by Overcamp and Fay, (1973). These 
authors pointed out that the concentrations of radiatively important trace species, such as water and 
ozone, could be greatly perturbed in aircraft wakes. In particular, a decrease in ozone would reduce 
the amount of UV absorption and direct heating of the plume, yielding a net cooling of the plume. 
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Similarly, enhanced concentrations of water or CO2 would increase the amount of IR cooling, and 
again yield a net cooling of the plume. Such net cooling could then lead to a net downward vertical 
velocity of the wake relative to the background atmosphere, as the perturbed air mass moves to 
maintain its potential temperature. The vertical velocity wouIcTbe given by 

w(km/day) = ~ (6) 

where w is the vertical velocity, Q (K/day) is the differential cooling rate of the wake relative to the 
background atmosphere, and G (K/km) is the so-called stability parameter, i.e., the vertical gradient 
in potential temperature. The value of G is fairly close to 10 K/km in the lower stratosphere. 

We have modified the radiative portion of our interactive 2-D model to calculate the expected 
differential cooling of aircraft wakes. We present calculations below of this differential cooling 
as a function of season and altitude, and for different assumptions about water enhancement and 
wake thickness. 

PHOTOCHEMICAL CALCULATIONS 

Calculations were carried out at an altitude of 18 km, for conditions representative of 47°N, 
June 1, starting at noontime. Our model of the wake regime will eventually be initialized from the 
results of the ARI model of the vortex regime. Since the ARI team has concentrated in modeling 
the jet regime, we have taken their calculated values for NO, N0 2 , NO3, N 2 0 5 and HN0 3 at a 
distance of 2000 ft from the nozzle, and simply diluted them to obtain a factor of 100 dilution from 
the nozzle to the beginning of the wake regime. Values for ozone, HO x species, as well as other 
calculated species were initialized to those calculated in our 2-D model for the above latitude and 
season. No dilution was applied to ozone, since the results of the ARI calculations at 2000 ft are 
very similar to the ozone concentrations in the background atmosphere. We also did not dilute OH 
and H0 2 , since these species have very short photochemical lifetimes. Adopted initial conditions 
for NO + N0 2 , HNO3, N 2 0 5 , and 0 3 are given in Table 1, as well as values for the background 
atmosphere taken from our 2-D model. 

1. The first hour 

a. Description of calculations 

The chemical evolution of different trace species was first calculated for a period of 1 
hour, with a time resolution of 50 seconds. This period of time was used to explore the sensitivity 
to heterogeneous processes, particularly on contrails, since it is unlikely that ice particles would 
last longer than one hour. 

Entrainment of background air into the wake air mass can occur through turbulent 
diffusion controlled both by dissipation of the internal kinetic energy and the turbulence of the 
background atmosphere. There is very little information on these processes for conditions in the 
lower stratosphere. The calculations for the first hour assumed a constant value for the entrainment 
parameter: 

dinVjt) 1 

dt t ' {) 

We have adopted a value of 780 seconds for t, which corresponds to an expansion of the volume by 
a factor of 100 in 1 hour. Calculations in the next subsection show results for slower entrainment 
rates. 
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The cases considered in our calculations are summarized in Table 2. Cases 1 and 2 
assume only gas-phase chemistry. Although large amounts of water are emitted by the aircraft, 
condensation into contrails and subsidence could remove a large portion of this water. We have 
therefore considered two extreme possibilites: all the water emitted by the engine s been 
condensed out, and the water mixing ratio is equal to that of the background atmosphere 'ase 1), 
or the water emitted by the engine stays in the wake, and the mixing ratio is given by the ? opriate 
dilution factor (Case 2). 

The heterogeneous reactions (2) - (5) are included in Cases 3 and 4 at rate similar to 
those utilized in our modeling of the Antarctic stratosphere (Rodriguez et a L, 1989). ’ . Case 3 we 
assume that no contrails are formed, and heterogeneous chemistry occurs only in the background 
sulfate aerosol layer presumably already entrained into the wake. In Case 4 we assume that all the 
water emitted by the aircraft has been condensed into cloud particles of 1 /im radius. We assume 
an initial concentration of 10 5 cm -3 for concentration nuclei, of which 1% are activated. With 
these assumptions, the available reacting surface concentration in the contrails is about a factor of 
10 4 higher than in the background sulfate aerosol layer, and heterogeneous reactions could occur 
very fast. We assume that contrails last for the duration of the calculations, but allow the density 
of particles to decrease with time in inverse ratio to the increase in plume volume. 

The reaction of N 2 Os with water (reaction 4) could be an important contributor to the 
conversion of NO x to HNO3. The rate of this reaction is however limited by the rate of conversion 
of NO2 to NO3 through the reaction 

N 0 2 + O3 — + NO3 + 0 2 , (8) 

which is very slow for stratospheric temperatures. Experiments of uptake of N0 2 on liquid surfaces 
(Lee and Tang, 1988) indicate an uptake coefficient of 6x 10 -4 , although this value could be even 
higher. If uptake of N0 2 could lead to reaction (8) occurring heterogeneously, then the conversion 
of N0 2 to NO3, N 2 0 5 and eventually HNO3 could be accelerated. Cases 5 and 6 were chosen 
to study this possibility. We assume a reaction efficiency of 6x 10~ 4 for reaction (8) at an ozone 
mixing ratio of 2.9 ppmv. Furthermore, we assume that this reaction efficiency is proportional to 
the ozone mixing ratio. In Case 5, the reaction is assumed to take place on bake ground sulfate 
aerosols, while in Case 6 the reaction takes place on contrails. 

b) Results 

The first question we seek to address is how much repartitioning of nitrogen species 
occurs in the wake regime due to photochemistry. The results in Table 3 shed light on this issue. The 
table shows the ratio of calculated NO x (NO x = NO + N0 2 + NO3 + 2xN 2 0 5 +H0 2 N0 2 +CINO3) 
and HNO3 after 1 hour to that calculated if no chemistry took place, and mixing ratios of the above 
trace gases were only affected by entrainment of background air. We see that for the first three 
cases the results are essentially the same. Therefore, there is little photochemical repartitioning 
of nitrogen species if we adopt only gas-phase chemistry, or even with heterogeneous chemistry 
occurring only on background sulfate aerosols. On the other hand, considerable repartit: rning 
occurs if we assume heterogeneous chemistry occurring on contrails (Cases 4, 5, and 6), with final 
concentrations of HNO3 enhanced by up to a factor of 2 over the entrainment-only values. 

The calculated time development of NO x and HNO3 is shown in Figures 1 and 2, 
respectively, for Cases 4, 5 and 6. The calculated concentrations of OH, NO x , HNO3, and O3 at 
the end of one hour are compared to values in the background atmosphere in Table 4. Inclusion of 
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heterogeneous chemistry on contrails leads to an enhancement of HN0 3 over background values, 
even if we allow for entrainment. Furthermore, heterogeneous conversion of N0 2 to N0 3 leads to 
further enhancement of HN0 3 and reductions in NO x during the first fifteen minutes; if this reaction 
takes place on contrails, almost compete conversion occurs in the first few minutes (Case 6). We 
also note that, except for Case 6, NO x concentrations are still enhanced relative to the background 
after 1 hour of the wake regime. This is due to the large initial concentrations of NO x ; further 
entrainment and mixing in the mesoscale regime will eventually bring the NO x concentrations to 
background values. 

An interesting feature of the photochemistry in the wake regime is the behavior of OH. 
For Cases 1 - 3, the large concentrations of NO x lead to reductions in OH relative to the background, 
due to the reaction of OH with N0 2 . In the presence of heterogeneous chemistry on contrails, 
however, the OH concentrations are enhanced. This is due to an increase in the HO x source due 
primarily to reactions (2) and (4), followed by photolysis of HOC1 and nitric acid. Such enhanced 
concentrations of OH accelerate the conversion of NO x to HN0 3 . In Case 6, rapid conversion of 
N 2 0 5 and C1N0 3 to HN0 3 reduces the magnitude of this source and the OH enhancement. 

The calculated time development of ozone in th e parc el is illustrated in Figure 3 for 
Cases 4-6. Ozone decreases during the first minute due to enhanced NO x and IIO x cycles. The 
decrease in ozone concentration is at most 15% during this period of time. For the rest of the hour, 
entrainment from the background atmosphere bring the ozone values to background levels by the 
end of the hour. Direct removal of ozone in the wake will not contribute to ozone reductions on a 
global scale. 

2. Scale-dependent Dispersion: The first two days 
a. Description of Calculations 

Our previous calculations assumed a constant, scale-independent entrainment time 
constant. We have also implemented a scale-dependent dispersion, following the formalism of 
Gelinas and Walton (1974). In this formalism, the wake is treated as a Gaussian plume, and its 
time-dependent volume, F(f), is assumed to expand according to the expression: 

dinV(t) K h , K v 

1 , 9*1 5 


In (9), the vertical and horizontal dimensions a v , ah are given by the equations: 


d 4= 2K ^ 

4 = 2 * •« 


( 10 ) 


The horizontal and vertical eddy coefficients A'*, I\ v can be assumed constant, or be related to 
the dimensions by expressions consistent with similarity theory. In particular, Gelinas and Walton 
suggest the expression: 

11 


K = 


be 1 / 3 a 4 / 3 I\ 


-l 


( 11 ) 
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where b is a constant of order 1, e the energy dissipation rate, and K m an asymptotic value for the 
eddy coefficient appropriate for large wake sizes. 

It is difficult to obtain reliable estimates for the parameters in equation (11). The 
large-scale distribution of trace gases in the stratosphere can be obtained with horizontal “eddy” 
coefficients of order 10 9 - 10 10 cm 2 s -1 (Ko et ai., 1985). At the same time, the length scale of the 
eddies responsible for these processes is much larger than that acting on aircraft wakes. Arguments 
from similarity theory suggest that horizontal eddy coefficients should be about as much as factor 
of 10 4 smaller for scales less than 1 km (Overcamp and Fay, 1973). Vertical eddy coefficients in 
the lower stratosphere are of order 10 3 cm 2 s -1 , and are expected to be maintained by small-scale 
processes. There is virtually no information on the energy dissipation parameter e. 

We carried out calculations of the time-development of aircraft wakes for a period 
of two days, using the initial conditions and chemistry of Case 2 (ie., enhanced water, but no 
heterogeneous chemistry), and with the following assumptions for the entrainment: 

i. No entrainment 

ii. “Slow” Entrainment: K ym - 10 4 cm 2 s _1 , e = 10~ 2 ergs/sec; K zm = 10 3 cm 2 s -1 . 

iii. “Fast” Entrainment: As above, but with K ym = 10 8 cm 2 s -1 .Our choices for 
the “fast” and “slow” entrainment are similar to those in the Gelinas and Walton 
formalism. Since heterogeneous chemistry on contrails is not expected to occur over 
several days, we did not include Heterogeneous processes in this set of calculations, 
but they can be easily incorporated in future calculations. 

b. Results 

The calculated entrainment time constant t (Equation 7) is shown in Figure 4 as a 
function of time, for the assumptions of slow (4a) and fast (4b) entrainment. Time constants for 
entrainment in the “slow” case increase from about 1 day at the beginning of the wake period 
to about 4 days after the 2-day calculation. The initial values for the time constant are about 30 
minutes in the “fast” entrainment case, closer to the assumptions in our previous set of calculations. 
However, given the faster increase in wake size in this case, the final values for the time constants 
after two days are similar to the “slow” case. It must be emphasized that this agreement is a 
coincidence due to our choice of parameters and integration times. 

Our calculations obtain dilution factors of 3 for the “slow” entrainment calculation, and 
20 for the “fast” entrainment calculation after two days. Given the initial concentrations adopted 
in Table 1, we see that significantly perturbed conditions could still exist two days after the passing 
of an aircraft, with our choice of dispersion parameters. This is further illustrated in Figures 5-8, 
which show the calculated time development of HNO3 (Figure 5), NO x (Figure 6), ozone (Figure 
7) and water (Figure 8) for cases with no entrainment, and “slow” and “fast” entrainment. 

The calculated HNO3 and NO x (Figures 5 and 6) are within percents of what would 
be calculated if only entrainment were acting on these species. The conclusion from the two-day 
calculation is therefore similar to the one over the first hour: in the absence of heterogeneous 
chemistry, there is little chemical conversion of NO x to HNO3 in the wake. Concentrations of NO x 
are still enhanced relative to the background, even after two days in the fast-entrainment case. 

The maximum calculated ozone reduction is slightly larger in this set of calculations 
(25%) than in the 1-hour results, since our entrainment time constants are slower. Such reductions 
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can be maintained for 1 day or longer if “slow” entrainment takes place. The ozone quickly 
recovers to background values in the “fast” entrainment case after two days. The “no entrainment” 
calculations show ozone still decreasing, albeit at a slower rate than during the first 100 seconds. 
The slowing down of ozone loss is due to the conversion of ozone-removing NO to other forms of 
NO x , such as NO2 and N2O5. 

Water concentrations are assumed to be initially about a factor of 70 higher than 
background. Figure 7 shows that enhancements of water of between factors of 5 and 20 over the 
background value (about 4 ppmv) could be still observed 2 days after the passage of the aircraft. 
The results in Figures 6 and 7 thus indicate that large perturbations in radiatively important species 
could be maintained for several days. These conditions could have important consequences for the 
net displacement and final deposition of aircraft effluents, as discussed in the next section. 

RADIATIVE CALCULATIONS 

We have adapted the radiative subroutines in our 2-D interactive model (Schneider et al., 
1989) in order to assess the differential cooling of aircraft wakes due perturbations in radiatively 
active species. This code includes a parameterization of solar heating from absorption in the UV 
and near infrared by ozone and water, as well as cooling for CO2, O3 and H 2 0 in the infrared. 
Parameterizations for overlapping bands are described in Wang and Ryan (1983). Calculated total 
heating rates with this code compare favorably with other calculations (Dopplick, 1970), although 
it is difficult to directly compare the contributions by different bands of species due to differences in 
adopted water and ozone profiles, and the apportioning of heating and cooling to different species 
in the spectral regions where bands overlap. 

Our calculations take the background atmospheric densities and UV and IR fluxes from our 
interactive 2-D model, since we expect that the aircraft wake will not alter the large-scale radiation 
field in the atmosphere. We assume that the wake temperature is the same as that of the background 
atmosphere. This may not be the case if dispersion of initial plume heat has not achieved completion; 
as a matter of fact, a relationship between the wake water vapor concentration and temperature can 
be obtained (Miake-Lye et al., 1991), indicating that temperatures at the beginning of the wake 
regime could be a few degrees higher than that of the surrounding atmosphere. We will include in 
our future work assessment of the effect of temperature differences on calculated cooling rates. 

We have carried out calculations of the differential cooling rates at altitudes between 13 and 
25 km, and for all latitudes and seasons. Modeling of the jet and wake regime and estimates 
of dilution factors indicate that concentrations of water could be as high as 300 ppmv (a factor 
of 70 enhancement over the background) at the beginning of the wake regime, in the absence of 
condensation in contrails (See Table 1 and discussion above). The results of our photochemical 
calculations indicate ozone reductions ranging from 15% to 25% (Figures 3 and 7). Based on 
this results, we have assumed for simplicity in our calculations water enhancements ranging from 
factors of 70 to 10, and ozone reductions of 20% over backgound values (2.9 ppmv). Although 
our code also includes the effects of C0 2 , preliminary estimates indicate that the impact of its 
enhancement would be minor compared to those of water and ozone. 

Another parameter of importance is the thickness of the wake. Our original code was adapted for 
incorporation into our 2-D model, with a vertical altitude resolution of about 3.5 km. Calculations 
assuming an enhancement of 70 in water over 3.5 km indicate that some water bands saturate, and 


54 


a “cooling to space” approximation is not appropriate. Since such vertical extents do not seem 
realistic for a single wake (although such vertical extent could be typical of a flight corridor), we 
further modified our code to allow for insertion of a wake of arbitrary vertical thickness at a given 
level. 

Our first set of calculations assumed a wake of 3.5 km thickness, with an enhancment in water 
by a factor of 70 and 20% reduction in ozone. This case is not considered to be realistic, and 
represents an upper bound to the expected enhancement in water over a flight corridor. Results for 
differential cooling rates (K/day) are presented in F igures 9-12 for days 30 (Figure 9), 120 (Figure 
10), 210 (Figure 11), and 300 (Figure 12). As discussed above, values for the vertical velocities 
can be obtained by dividing the numbers in Figures 9-12 by a factor of 10 (Equation 7). Vertical 
velocities as large as 0.4-0.5 km/day could be obtained for wakes deposited at altitudes above 22 
km. Even at lower altitudes, vertical velocities of 0.2 km/day are obtained at mid-latitudes. 

Calculations in Figures 1 3- 16 are carried out for the same days and the same water enhancement 
as in the previous set, but assuming a vertical thickness of 500 m for the wake. Comparison of 
these two sets of calculations illustrates the effect of water band saturation. These effects are less 
pronounced for a wake of only 500 m thickness, and thus the calculated cooling rates (and vertical 
velocities) are actually larger. Vertical velocites close to 1 km/day are obtained at altitudes above 
22 km, with about 0.5 km/day being calculated at lower altitudes. 

The gas-phase concentrations of water assumed in the above calculations may be difficult to 
maintain without any condensation (see e.g., Miake-Lye et a i., 1991). We have thus carried out a 
third set of calculations, assuming only a factor of 10 enhancement in water. Results are shown in 
Figures 17-20. Vertical velocities ranging from 0.2 to 0.3 km/day are obtained in this case. We note 
that, even for a 500 m wake, the IR cooling rates do not scale linearly with water concentration, 
suggesting that saturation effects are still operant over this range. Work is in progress to further 
explore this issue. 

The results presented in Figures 9-20, when combined with the time-development of water in 
Figure 8, indicate that the aircraft wake could sink 1-2 km during the first two days. These results, 
however, depend on several processes and parameters which are at present not well quantified. 
These include: 

i. Formation and lifetime of contrails, as well as possible removal of water by contrail 
subsidence. 

ii. Time constants for turbulent dispersion of the wake. 

iii. Other processes contributing to wake dispersion and deformation. In particular, the impact 
of wind shear in reducing the smallest dimensions of the wake to a size when molecular 
diffusion becomes important should be considered (Prather and Jaffee, 1990). 

DISCUSSION AND CONCLUSIONS 

We have developed a Langrangian box model to study time development of trace species in 
the wake regime of an aircraft plume. This model has been used in preliminary calculations 
in order to assess the amount of conversion of NO* to HNO3. Our calculations corroborate 
the CIAP results of the mid-70’s (Gelinas and Walton, 1974; Hoshizaki et al., 1975) i.e., that 
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there is very little conversion of NO x to HNO 3 if only gas-phase chemistry is assumed to occur. 
Substantial conversion could occur, however, if heterogeneous chemistry takes place on contrails 
which could provide large enhancement in available surfaces for heterogeneous reactions. In 
particular, heterogeneous conversion of NO 2 to N0 3 could accelerate the NO x to HNO 3 conversion 
by heterogeneous reactions. Heterogeneous chemistry, particularly reactions (2) and (4), provide 
an additional source of HO x and lead to enhanced concentrations of OH in the wake. 

We stress that these calculations presented here are of an illustrative nature. The duration of 
heterogeneous chemistry has been exaggerated, since we do not expect contrails to persist for one 
hour. There are also varying degrees of uncertainty in the rates of the heterogeneous reactions (2) - 
(4), depending on the composition of the reacting surface. The heterogeneous conversion of NO 2 
to HNO 3 has not been experimentally established. 

We have also extended our calculations to a time-period of two days, and investigated the sen- 
sitivity of our results to different assumptions about scale-dependent dispersion. Our calculations 
indicate that a chemically-perturbed environment could persist for two days after passage of the 
aircraft. In particular, concentrations of radiatively important species such as ozone and water 
could be reduced by about 20% and enhanced by factors of 4-70 respectively during this period of 
time. 

The above results motivated an examination of the radiative balance in aircraft wakes. Enhance- 
ments in water and reductions in ozone yield differential cooling rates relative to the background 
which could be as large as 10 K/day. Calculated downward vertical velocities range from 0.2 to 1 
km/day, depending on the altitude of wake deposition and the assumed water concentration. The 
net vertical displacement of the wake will be a sensitive function of the rate of dispersion of the 
wake. 

We stress that the implications of both our photochemical and radiative results depend crucially 
on the adopted rates of dispersion. Information about this parameter is at present practically non- 
existent. Furthermore, other processes not directly included in this formation could be equally 
important. Prather and Jaffee (1990) examined the dispersion of chemically-perturbed polar air 
masses due to random shear in a turbulent field. This model actually predicts a decrease with 
time in the smallest dimension of the air mass, leading to an elongated “ribbon” which probably 
wraps around itself. When the smallest dimension of this ribbon reaches values of tens of meters, 
molecular diffusion becomes important, and the perturbed air mass mixes with the background in 
a matter of hours. 

The above picture presents an alternative and probably more realistic view of wake dispersion. 
The “volume” of the wake corresponding to equation (1) could be considered to be that containing 
all of the “ribbon” wrapped around itself. At the same time, this model allows for patchiness 
inside the wake, and thus dispenses with the assuption of a Gaussian profile. This wake patchiness 
could presumably alter the photochemical calculations. Modeling of this kind of process requires 
a different approach from the one utilized in this and previous studies, and could lead to faster time 
constants than assumed here. Still, it must be noted that the rate of random strain by a turbulent 
stratospheric wind field is still an unknown parameter; estimated time constants range from 1 - 10 
days (Prather and Jaffee, 1990). " 

Our calculations point to areas which deserve further study. These include: 
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a. Are contrails formed in the wake regime, and what is the expected lifetime? 

b. Do the enhanced concentrations of HN0 3 increase the probability of forming nitric acid 
tri-hydrate contrails? 

c. If contrails are formed, how big are the particles, and could their subsidence contribute to 
a redistribution of HNO 3 and water? 

d. What are the entrainement rates in the lower stratosphere, and how would they compete 
with heterogeneous chemistry and contrail formation? 

e. What are the roles of random strain and molecular diffusion in the final mixing of aircraft 
wakes, and what are the time constants for these processes? 

REFERENCES 

Dopplick, T. G. (1970) Global radiative heating of the Earth’s atmosphere. Planetary Circulation 
Project Report 24. Department of Meteorology, MIT, Cambridge, MA. 

Gelinas, R.J., and J.J. Walton (1974) Dynamic-kinetic evolution of a single plume of interacting 
species. J. Atmos. Sci., 31, 1807. 

Hoshizaki, H., L.B. Anderson, R.J. Conti, N. Farlow, J.W. Meyer, T. Overcamp, K.O. Redler, 
and V. Watson (1975) Aircraft wake microscale phenomena. Chapter 2 in CIAP Monograph 3, 
Department of Transportation, Washington, DC, DOT-TST-75- 53. Ko, M. K. W., K. K. Tung, 
D. K. Weisenstein, and N. D. Sze (1985) A zonal mean model of stratospheric tracer transport in 
isentropic coordinates: Numerical simulations for nitrous oxide and nitric acid. J. Geophys. Res., 
90, 2313. 

Lee, J.H. and I.N. Tang (1988) Accomodation coefficient of gaseous NO 2 on water surfaces. 
Atmos. Environ. 22, 1147-1151. 

Leu, M.-T. (1988a) Laboratory studies of sticking coefficients and heterogeneous reactions 
important in the Antarctic stratosphere. Geophys. Res. Lett., 15, 17. 

Leu, M.-T. (1988b) Heterogeneous reactions of N 2 O with H 2 O and HC1 on ice surfaces: 
Implications for Antarctic ozone depletion. Geophys Res. Lett., 15, 851. 

McElroy, M.B., R.J. Salawitch, S.C. Wofsy, and J.A. Logan (1986) Antarctic ozone: Reductions 
due to synergistic interactions of chlorine and bromine. Nature, 321, 759. 

Miake-Lye, R.C., M. Martinez-Sanchez, R.C. Brown and C.E. Kolb (1992) Plume and wake 
dynamics, mixing, and chemistry behind an HSCT aircraft. J. Aircraft, accepted for publication. 

Molina, M.J., T.L. Tso, L.T. Molina and F.C.Y. Wang (1987) Antarctic stratospheric chemistry 
of chlorine nitrate, hydrogen chloride and ice: Release of active chlorine. Science, 238, 1,253. 

Mozurkewich, M. and J.G. Calvert (1988) Reaction probability of N 2 0 5 on aqueous aerosols. 
J. Geophys. Res., 93, 15,889. 

Overcamp, T. J., and J. A. Fay (1973) Dispersion and subsidence of the exhaust of a supersonic 
transport in the stratosphere. J. Aircraft, 10, 720-728. 

Prather, M.J., and A.H. Jaffe (1990) Global impact of the antarctic ozone hole: Chemical 
propagation. J. Geophys. Res., 95, 3473- 3492. 


57 


Schneider, H. R., Malcolm K. W. Ko, N. D. Sze, G.-Y. Shi, and W.-C. Wang (1989) An 
evaluation of the role of eddy diffusion in stratospheric interactive two-dimensional models. J. 
Atm. Sci., 46, 2079-2093. 

Rodriguez, J.M. M.K.W. Ko and N.D. Sze (1988) Antarctic chlorine chemistry: Possible global 
implications. Geophys. Res. Lett., 15,257. 

Rodriguez, J.M. , M.K.W. Ko, and N.D. Sze, S.D. Pierce, J.G. Anderson, D.W. Fahey, K. 
Kelly, C.B. Fanner, G.C. Toon, M.T. Coffey, L.E. Heidt, W.G. Mankin, K.R. Chan, W.L. Starr, 
J.F. Vedder, and M.P. McCormick (1989) Nitrogen and chlorine species in the spring Antarctic 
stratosphere: Comparison of models with Airborne Antarctic Ozone Experiment observations. J. 
Geophys. Res., 94, 16,683. 

Solomon, S„ R.R. Garcia, F.S. Rowland, and D.J. Wuebbles (1986) On the depletion of antarctic 
ozone. Nature, 321, 755. 

Tolbert, M.A., MJ. Rossi, R. Malhotra, and D.M. Golden (1987) Reaction of chlorine nitrate 
with hydrogen chloride and water at anarctic stratospheric temperatures. Science, 238, 1258. 

Tolbert, M.A., M.J. Rossi and D.M. Golden (1988a) Heterogeneous chemistry related to Antarc- 
tic ozone depletion chemistry: Reactions of N 2 Os with H 2 0 and HC1 on ice surfaces. Science, 
240,1018. 

Tolbert, M.A., M.J. Rossi and D.M. Golden (1988b) Heterogeneous interactions of chlorine 
nitrate, hydrogen chloride, and nitric acid with sulfuric acid surfaces at stratospheric temperatures. 
Geophys. Res. Lett., 15, 847. 

Van Doren, J.M., L.R. Watson, P. Davidovits, D.R. Worsnop, M.S. Zahniser and C.E. Kolb 
(1991) Uptake of N 2 Os and HN0 3 by aqueous sulfuric acid droplets. J. Physical Chem. In press. 

Wang, W.-C., and P. B. Ryan (1983) Overlapping effect of atmospheric H 2 0, C0 2 and 0 3 in 
the C0 2 radiative effect. Tellus, 35B, 81-91. 

Wofsy, S.C. (1978) Temporal and latitudinal variations of trace gases: A critical comparison 
between theory and experiment. J. Geophys. Res., 83, 364-378. 


58 


Table 1: Adopted Initial Conditions 


Species 


NO + NO2 


Initial condition* 
(mixing ratio) 

4.48 x 10' 7 


HNO3 


2.80 x 10-3 
2.20 x 10-10 


Background atm. 


1.55 x 10-9 
4.75 x 10-9 
2.60 x 10-10 


N2O5 

O3 


2.91 x 10'S 


2.91 x 10-6 















Table 2: Cases Considered 









































Table 4: Comparison of Trace Species in Wake to Backgrou 

Concentrations (after 1 hour) 


Case 


1 


2 


3 


4 


5 


6 


OH 

OH(bckgd.) 

0.52 


0.63 


0.65 


2.83 


2.84 


0.80 


hno 3 

HNQ 3 (bckgd.) 

1.00 


1 .00 
1.07 
1 .39 
1.56 
2.30 


NO* 

NOx(bckgd-) 


03 

0 3 (bckgd.) 


2.61 


2.60 

2.60 


1.96 


1.69 


0.47 


0.998 

0.999 

0.999 


0.998 


0.998 


0.997 






































Mixing ratio 



Figure 2: Same as Figure 1, except for NO x . 
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Figure 4. Calculated entrainment time constant t (see equation 7) for the 
cases with scale-dependent diffusion as a function of time. Panel (a) 
shows the calculations for “slow” entrainment, while panel (b) denotes 
the calculated x for “fast” entrainment. See text for details. 





Figure 5. Calculated time development of nitric acid for the assumptions 
of case 2, except that a scale-dependent parameterization for entrainment 
is adopted. Results for no entrainment (solid line) and slow (dashed line) 
and fast (dotted line) scale-dependent dispersion are presented. 
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Mixing ratio 1.0000e-06 


Figure 



7. Same as Figure 5, except for calculated ozone. 
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Difference of net heating (for 30) 


°K/DAY 



Latitude 


Figure 9. Differential heating rate (°K/day) as a function of altitude and 
latitude, for day 30. We assume that the concentrations of water and 
ozone are respectively increased and decreased by factors of 70 and 0.8 
over the background values. These increases are assumed to hold over the 
whole grid vertical width (3.5 km in our 2-D model). 
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Figures 10. Same as Figure 9, except for day 120. 
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Figures 12. Same as Figure 9, except for day 300. 
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Figure 13. Same as Figure 9, except that the vertical thickness of the 
plume is taken to be only 500 m. 
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Total heating (case2— easel) day 210 
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Figure 15. Same as Figure 13, except for day 210. 
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Total heating (case2— easel) day 300 
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Figure 16. Same as Figure 13, except for day 300. 
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Total heating (case4— easel) day 30 

°K7DAY 



Figure 17. Same as Figure 13, except that the water is enhanced only by a 
factor of 10. 
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Total heating (case4— easel) day 120 
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Figure 18. Same as Figure 17, except for day 120. 
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Figure 20. Same as Figure 17, except for day 300. 







Chapter 4 

Laboratory Studies of Stratospheric Heterogeneous Chemistry 
INTRODUCTION 

The purpose of this aspect of our program is to evaluate the role of heterogeneous chemistry 
relating to the potential atmospheric effects of emissions from a fleet of supersonic aircraft flying in 
the stratosphere. Research needs in this area are motivated by a relative lack of knowledge of 
chemistry on or in aerosol particles compared to that of the homogeneous gas phase. Significant 
advances in our understanding of heterogeneous chemical mechanisms are required for quantitative 
prediction of stratospheric chemistry with or without added aircraft emissions. 

The importance of heterogeneous chemical reactions on stratospheric aerosol surfaces was 
originally demonstrated by measurements and models of chemical constituents in the 1987 
Antarctic spring (Brune et al., 1989; Fahey et al., 1989; Solomon, 1988; Salawitch et al., 1988). 
Subsequent arctic winter observations of CIO (Brune, 1989), NOy (Fahey, 1989) and polar 
stratospheric clouds (Hofmann et al., 1989b) reinforced the theory that heterogeneous reactions on 
PSCs are responsible for both denitrification and chlorine activation which are necessary 
precursors for the large ozone depletions observed in the Antarctic polar stratosphere. 

More recent analysis of Arctic data, directly correlating aerosol appearance and perturbation of 
(radical) chemical partitioning (Kawa et al., 1991a, 1991b), serves to highlight the complexity of 
the complete gas/aerosol system. Those studies have focused on discrepancies between Arctic and 
Antarctic aerosol observations and the H2O/HNO3 phase diagram (Hanson and Mauersberger, 
1988a) that is believed to control aerosol appearance. These observations of apparent chemical 
perturbation via aerosols have been supported by a growing collection of laboratory data which 
have shown that heterogeneous processes are more rapid than previously believed and that 
chemical conversions on cold aqueous surfaces can result in gaseous products (Molina et al., 1987; 
Tolbert et al., 1987, 1988a, 1988b; Quinlan et al., 1990; Reils et al., 1990; Worsnop et al., 1988; 
Watson et al., 1990; Van Doren et al., 1990, 1991; Leu, 1988a, 1988b; Moore et al., 1990; 
Mozurkewich and Calvert, 1988). These data, coupled with the field measurements in the Arctic 
and Antarctic, provide a qualitative picture of the role of heterogeneous chemistry in the polar 


stratosphere. 

Through the following reactions occurring on PSCs 

HCl(g) + H 2 0(s) ~> HCl(s) R1 

HN0 3 (g) + H 2 0(s) --> HN0 3 (s) R2 

N 2 0 5 (g) + H 2 0(s) ~> 2HN0 3 (s) R3 

ClN0 3 (g) + H 2 0(s) ~> HN0 3 (s) + HOCl(g) R4 

ClN0 3 (g) + HCl(s) ~> HN0 3 (s) + Cl 2 (g) R5 

N 2 0 5 (g) + HCl(s) -> HN0 3 (s) + ClN0 2 (g) R6 


the chlorine reservoir species C1N0 3 and HC1 are converted to gas phase products which photolyze 
to eventually form ozone-reactive CIO. Denitrification of the stratosphere also occurs from the 
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conversion of CINO3 and N2O5 (and condensation of gaseous HNO3) to HNC>3(s) which impedes 
the formation of CINO3, thereby maintaining high levels of CIO. 

A key question is the formation and chemical composition of the aerosol particles referred to as 
(s) in R1-R6. It was the presence of PSCs that correlated with the spring ozone depletion over 
Antarctica. At the coldest temperatures observed there (185 K), PSCs clearly are composed of ice 
particles, probably containing significant amounts of HNO3. At least some of these are large 
enough to precipitate, contributing to extreme denitrification and dehydration of the Antarctic 
stratosphere. Enhanced aerosol densities were also observed at higher temperatures, above the 
frost point of ice (McCormick et al., 1982). These were predicted to be composed of nitric acid 
trihydrate (NAT) solid solutions (Toon et al., 1986). Comparison of stratospheric mixing ratios of 
HNO3 and H2O with laboratory measurements of temperature dependent partial vapor pressures 
above H2O/HNO3 mixtures were consistent with this (Hanson and Mauersberger, 1988a, 1988b). 

The results presented in this report have confirmed and extended the Hanson and Mauersberger 
(1988a, 1988b) results. In addition to reproducing their equilibrium vapor pressure diagram for 
NAT, our results have resolved non-equilibrium effects relating to slow diffusion within the 
condensed solid phase and to supersaturation of gas phase HNO3 with respect to NAT nucleation. 
In the stratosphere, slow solid diffusion may result in non-homogeneous NAT/ice particles, 
especially in air parcels that cycle through cooling and warming events (Wofsy et al., 1990). As 
mentioned above, HNO3 supersaturation can explain discrepancies between Arctic observation of 
aerosols as predicted by the equilibrium NAT phase diagram (Kawa et al., 1991b). 

Over the past few years, these field and laboratory observations have spurred attempts to 
model the potential role of heterogeneous chemistry on a global scale in the stratosphere. Hofmann 
and Solomon (1989) estimated that the influx of sulfate aerosol into the stratosphere from the El 
Chichon volcanic eruption in 1983 may explain reduced NO2 and O3 levels observed in the 
northern hemisphere in the subsequent two years. Others have modeled potentially significant 
effects on global ozone levels by extrapolating simple parameterizations of (polar) heterogeneous 
chemistry to mid-latitudes (Rodriguez et al., 1988; Isaksen, 1989). At the same time, the diversity 
of in situ chemical and aerosol observations both in general and especially during the arctic winter 
have served to emphasize the complexity of the phenomena. 

The problem with these modeling attempts is that they rely on rather crude parameterizations of 
heterogeneous reaction rate parameters based on the limited laboratory data available for reactions 
R1-R6. For example, Hofmann and Solomon (1988) and Isaksen (1989) used values from 
previous work in our laboratory for N2O5 reaction (R3) and reactivity coefficient data of Tolbert et 
al. (1988b) for CINO3 with and without HC1 (R4 and R5) on sulfuric acid surfaces. Rodriguez et 
al. (1988) parameterized the effective reaction rates for R1 and R6 in terms of the equivalent gas 
phase reaction rate necessary to explain the high CIO levels observed in the Antarctic spring and 
then scaled those rates to the global stratosphere by scaling according to the total aerosol surface 
area. 

Such parameterizations ignore both the diversity of the aerosol conditions in the stratosphere 
and the details of the experiments that measured the heterogeneous rates. Each experiment 
measures an effective heterogeneous rate that can depend on parameters such as trace gas density, 
gas-surface exposure time, and method of surface preparation. In our laboratory kinetic 
measurements we have shown that experimental rate determinations are often only a measure of 
gas solubility in the condensed phase. Likewise, our new vapor pressure measurements indicate 
that surface composition can deviate significantly from the interior solid, depending on gas/surface 
interactions (Worsnop et al., 1992). For reactive processes such as R5-R6, detailed knowledge of 
the chemical mechanism is needed to quantitatively relate laboratory kinetic data to atmospheric 
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conditions. Furthermore, that modeling must take into account the distribution of aerosol density 
and composition that will in general be variable in H 2 O, HNO 3 and H 2 SO 4 . 

This need for more sophisticated treatment of both the atmospheric surfaces and the 
heterogeneous chemical kinetics is the backdrop on which potential heterogenous pertur' tions due 
to stratospheric aircraft need to be considered. 


PERTURBATIONS FROM AIRCRAFT EMISSIONS 

The HSRP program is indeed focusing much of its effort on development of stratospheric 
modeling programs, since the effect of aircraft emissions can only be predicted via model 
simulation. Clearly the first step involves improving the capability of modeling the present 
stratosphere. Along those lines, appropriate inclusion of heterogeneous chemistry is a major issue 
for model development In fact, this has been demonstrated in a new calculation by the group at 
AER (Weinsenstein, et al., 1991), which has included reaction R3 for N 2 O 5 on global background 
H 2 SO 4 aerosol. Using reaction rates for R3 measured in our laboratory (Van Doren et al., 1991) 
and others (Mozurkowich and Calvert, 1 988), they predict that R3 shifts partitioning of NO x to 
HNO 3 , which in turn shifts Cl x /NO x chemistry leading to O 3 depletion by Cl x . This calculation 
actually predicts that NO x emission from HSCT aircraft would increase O 3 by partially neutralizing 
the Cl x destruction cycle. At this point, such a calculation re-emphasizes the need to treat 
stratospheric heterogeneous chemistry as completely as possible. 

The impact of emissions from stratospheric aircraft must be considered at local and global 
levels. The aircraft fleet is projected to fly at altitudes between 18 and 23 km. This corresponds to 
the lower stratosphere, extending down to the top of the tropopause and overlapping the natural 
sulfate aerosol layer. The local temperature and trace gas profiles of aircraft plumes can be 
expected to have spatial and temporal scales of hundreds of meters and hours. Emission of high 
(relative to background stratospheric levels) densities of H 2 O, NOx, SO 2 , and condensation nuclei 
can be expected to significantly perturb aerosol density and composition within these plumes. 
These plumes are expected to mix on a timescale of days, giving a "zonal average" in aircraft flight 
corridors. Obviously on a longer timescale there is dispersion throughout the stratosphere as well 
as mixing down into the top of the troposphere. 

Quantitative modeling of perturbed stratospheric profiles requires detailed coupling of both 
homogeneous and heterogeneous chemistry with transport models, which is the overall goal of the 
HSRP program. In this context, the plume mixing scenario represents an extension of the range 
of gas and temperature conditions that need to be investigated. In fact, the higher temperatures and 
trace gas densities associated with the aircraft plumes are somewhat easier to study in the 
laboratory than those of the unperturbed stratosphere. This may be particularly true for 
heterogeneous interactions of NO x and SO 2 with aerosols 

OVERVIEW OF EXPERIMENTAL PROGRAM 

In separate collaborations with Boston College and Harvard University, Aerodyne Research 
has developed laboratory projects for measuring heterogeneous kinetic and thermodynamic 
parameters. Together with the group of Professor Paul Davidovits of Boston College, a technique 
based on coupling a fast-moving train of liquid droplets with a gas flow tube has successfully 
measured a number of gas/surface accommodation/reaction processes of tropospheric and 
stratospheric interest. Together with the group of Professor Steven Wofsy of Harvard an 
apparatus consisting of a cooled infrared absorption cell has been used for the determination of 
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vapor pressures of solid solutions. Both of these projects provide state-of-the-art experimental 
approaches to the difficult problems associated with studying heterogeneous processes. 

At the same time, this collaborative approach provides a unusually broad perspective on the 
issues involved. Although the two experimental projects have been programmatically separate, 
they complement and synergistically interact on personnel, technical and scientific levels. The 
collaboration of Aerodyne staff members, Professors Davidovits and Wofsy, and students and 
postdoctoral associates from the academic groups has brought together a unique combination of 
physical and atmospheric scientists. On a technical level, the two projects utilize the same diode 
laser techniques for trace gas detection as well as bringing together different perspectives on the 
practical problems associated with preparing well-defined surfaces in the laboratory. Scientifically 
the issues of aerosol composition and chemistry obviously are deeply intertwined. 

HETEROGENEOUS KINETIC STUDIES 
Introduction 

Because of the small surface area and cold temperature of stratospheric aerosols, it was 
believed until recently that the likelihood of significant stratospheric heterogeneous chemistry was 
very small. Field measurements correlating Antarctic ozone depletion with PSCs and surprisingly 
large heterogenous rates measured in the laboratory have changed that view. However, a much 
more detailed understanding of the underlying chemical mechanisms is required before these 
results can be quantitatively applied to heterogeneous processes in the global stratosphere, with or 
without added perturbation from aircraft emissions. Below we discuss the issues involved in 
improving that understanding. 

The rate of uptake of gas into a condensed phase can be expressed as a pseudo first order rate 
by the expression 


ycA (l) 


where c is the molecular velocity (cm s-1) and A is the aerosol surface area density (cm2 cm-3). 

The parameter y, often referred to as the "sticking coefficient", is the probability of gas uptake into 
a surface upon collision with the surface. More precisely, in the case of direct accommodation, as 
for HC1 and HNO3 in R1 and R2, yis the mass accommodation coefficient; for uptake involving 
reaction, as for N2O5 and CINO3 in R3 and R4, yis the reactivity coefficient 

Analysis of most laboratory data and its application to atmospheric modeling has been solely in 
terms of Eq. (1). Experimental uptake measurements are simply compared to the surface collision 
rate to extract a value of y. This over-simplified interpretation ignores the fact that, even for simple 
accommodation processes such R 1 and R2, experimentally observed gas uptake measurements in 
general involve a combination of processes. In particular, one must always evaluate the degree of 
saturation of the surface. In the case of bimolecular processes such as R5 and R6, the dependence 
of heterogeneous rates on gas concentration is obviously important, especially when laboratory 
densities are much higher than in the atmosphere. 

Equation (1) does provide a simple framework to evaluate the potential role for heterogeneous 
processes. Under typical stratospheric conditions the surface area. A, peaks at about lO 8 cm2 
cm-3 between 15 and 20 km altitude, decreasing to about 10* 9 at 30 km (Turco et al., 1981). Thus, 
even in the limit of y = 1 , the maximum heterogeneous rate is on the order of 10 4 s* 1 . In other 
words, heterogenous processes in the stratosphere can be significant only for gaseous species with 
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long reactive lifetimes. Until recently it was generally believed that y would be small at 
stratospheric temperatures, making heterogeneous kinetics unimportant 

The polar field observations and various laboratory results have changed that view. The 
correlation between ozone depletion and the presence of PSC’s combined with laboratory 
measurements of y > 0.0 1 even for bimolecular heterogeneous reactions such as R5 and R6 now 
indicate that heterogeneous processes can indeed be significant, at least under polar conditions. 
More generally, results from our laboratory have shown that direct mass accommodation for 
molecules such as HC1 and HNO3 on aqueous surfaces (R1 and R2) is efficient (Watson et al., 
1990). The accommodation probability for a number of molecules has a negative temperature 
dependence so that y approaches unity at stratospheric temperatures. This extrapolation is 
consistent with results for H2O, HNO3 and HC1 accommodation on cold (200 K) ice and NAT 
surfaces which give lower limits on y approaching unity within experimental limits (Leu, 1988a; 
Moore et al., 1990; Molina et al., 1987; Tolbert et al., 1988; Quinlan et al., 1990). 

For such large values of y, the rates of direct accommodation processes such as R1 and R2 
become irrelevant. For small aerosol particles under stratospheric conditions, molecular density in 
the particles will equilibrate with the gas density on timescales much shorter than a second. Thus, 
for both atmospheric and laboratory conditions the important parameter is the gas solubility in the 
condensed phase. 

A case in point is the uptake of HC1 into aqueous sulfuric acid solution measured as a function 
of H2SO4 mole fraction. The decrease in observed HC1 uptake with increasing H2SO4 mole 
fraction may be modeled by the equation 

y = — (2) 


1 + — CLC — 

_Q_ 

1/2 

4 H* RT 

-t . 



The key parameters are the effective HC1 Henry’s law coefficient H* (which depends on the 
dissociation of HCl(aq) into H+ and C1-), the liquid diffusion coefficient D, the gas/liquid exposure 
time (2 x 10-3 s in our experiment), and the true mass accommodation coefficient a. The data can 
be quantitatively modeled using bulk solution parameters, consistent with uptake limited by 
time-dependent, diffusion-limited saturation of the liquid surface. This mechanism is also 
consistent with data of Reihs et al. (1990). Applying this analysis to their HC1 uptake data, one 
can extract a value of H*(D)-l/2 for HC1 in 65 wt% H2SO4 solution at 205 K. The key question 
there is the physical state of the aqueous H2SO4 surface, which determines the value of D. 
Depending on whether the surface is frozen or not, the value of H* for HC1 is 30 to 10 4 times 
larger at 205 K than at 283 K. 

The values of H* and D for HC1 are also critical to evaluating the rates of the bimolecular 
reactions R5 and R6. In analogy with Eq. (2), one can derive an expression for y that compares 
the rates of condensed phase reaction and re-evaporation due to saturation. In the limit of reaction 
rate limited uptake, 


7R = 4H ^ RT { Dk R H*RT [HCl(g)] ) 1/2 
c 


where Hr is the Henry ’s law solubility of reactant R and kR is the bimolecular rate constant for 
condensed phase reaction of HCI with reactant R. 
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The heterogeneous rates reported for reactions such as R5 and R6 (y> 0.01) need to be 
compared to Eq. (3). Again a key parameter is the diffusion coefficient for HCl(s). Originally D 
for HC1 in ice was reported to be unusually large, -10-5 cm 2 s- 1 (Molina, 1987), which indicated 
that reactions R5 and R6 could access all HCl(s) dissolved in aerosol particles. However, in a 
separate study it was pointed out that HC1 in single crystals of ice actually is slow, though HC1 is 
readily adsorbed and transported along grain boundaries (Wolff et al., 1989). Using a typical 
value of D = 10-U cm 2 s- 1 in solid solution, assuming H*(HC1) = 10 10 M atm-1, and assuming 
diffusion limited kR of 10 4 M* 1 s~ [ , for [HCl(g)j = 109 cm-3, we obtain 

y = Hr x 10‘6 

This estimation is based on upper limits for the bulk phase parameters. For example, recent 
measurements indicate that the actual equilibrium solubility of HC1 in ice may be significantly 
smaller (Hanson and Mauersberger, 1989). While little is known about the values of Hr for 
species such as CINO3 and N2O5, it is hard to reconcile the large observed rates for R5 and R6 
(y > 0.0 1). It is also significant that these rates appear to be more or less independent of the 
condensed surface, whether ice, NAT or aqueous H2SO4. This, combined with the fact that R5 
and R6 rates are comparable to those for R3 and R4 reactions with H20(s) that is in huge excess, 
points toward a surprisingly efficient surface reaction. 

Evaluation of surface reaction rates requires knowledge of surface concentrations of adsorbed 
reactant molecules. Turco et al. (1989) have discussed the thermodynamics of so-called adsorption 
isotherms and estimated that surface binding energies on the order of 10 to 20 kcal mol-1 are 
required to give significant surface densities. In their analysis, they implicitly concluded that the 
reported efficiencies for reactions R3 to R6 were consistent with such surface densities. The rate 
of reaction of two adsorbed species such as HC1 and CINO3 also is dependent on their mobility on 
the surface, which is related to the entropy change associated with the adsorption isotherm. It 
should be noted that in general the rate of such a bimolecular surface reaction is proportional to the 
product of the molecular surface densities, each of which is in general proportional to the 
respective gas densities. However, in the limit of high (e.g. laboratory) gas densities it is easy to 
saturate surface densities and observe heterogeneous rates that are independent of gas density. 

This mechanism of surface adsorption is consistent with results from our laboratory for mass 
accommodation on aqueous surfaces. For the series of molecules H2O2, HNO3, HC1 and N2O5, a 
negative temperature dependence was observed that was consistent with surface binding energies 
of at least 5 to 10 kcal moH (Worsnop et al., 1989; Van Doren et al., 1990). In the case of SO2, 
the density of the surface state was directly measured, again consistent with a surface binding 
energy on the order of 10 kcal moH (Jayne et al., 1990). These results were interpreted in terms 
of general accommodation mechanism involving initial adsorption on the liquid surface. The 
rate-limiting step for accommodation into the liquid was then the rate of solvation of the 
surface-bound molecule. The actual accommodation coefficient then reflected the competition 
between solvation and re-evaporation of the surface state. 

These measurements need to be extended to lower temperatures for the chemical species and 
surfaces relevant to R1 - R6 in the stratosphere. If indeed heterogeneous chemistry occurs via 
surface adsorption, the surface densities will be observable using the experimental techniques 
developed in our laboratory. The key capability is the millisecond resolution for gas/liquid 
exposure times. This permits the resolution of the initial adsorption from the equilibration of 
adsorption on longer timescales. 

In this report we present experimental data on the uptake of HC1 on cold sulfuric acid droplets. 
Understanding the uptake, coefficients in terms of the solubility of HC1 is the first step in probing 
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surface reactivity of gas phase species on stratospheric sulfate aerosols reased by adsorption on 
aerosols. 

Experimental Measurements of HC1 Uptake On Sulfuric Acid 
Droplets 

We have measured the uptake of gas phase HC1 by droplets of sulfuric acid as a function of 
droplet temperature and concentration. A fast-moving (~3500 cm/s) train of concentrated sulfuric 
acid droplets (40, 50, and 70 wt%) is formed by passing pressurized acid through a 70 pm orifice 
in a thin platinum disk that is mechanically vibrated using a piezoelectric ceramic. The droplet 
surface area is varied by switching the droplet frequency which changes the diameter and spacing 
of the droplets (droplet diameter = 160 pm). The droplet train travels coaxially down a 30 cm 
long, 1.5 cm ID glass flow tube in which helium is used as the carrier gas. HC1 mixed with 
helium is introduced using one of three Teflon loop injectors positioned downstream of the orifice. 

The contact distances of HC1 with the droplet stream are 5.5, 16.0, and 29.5 cm which, for a 
droplet velocity of 3500 cm/s, correspond to droplet-HCl exposure times of 1.6 - 8.4 ms. A White 
cell connected to the downstream end of the flow tube allows for tunable diode laser infrared 
absorption measurements of the HC1 concentration as a function of the total surface area of the 
stream exposed to HC1. 


Temperature Controlled 



Figure 4.1. Schematic of the Apparatus for the Measurement of Uptake Coefficients 
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The acid is cooled to temperatures in the range 230 - 260 K using a refrigeration bath and a 
circulation pump. A thermocouple placed just above the nozzle is used to monitor the liquid 
temperature. Tunable diode laser infrared absorption measurements of the ambient water vapor in 
the flow tube are used to determine droplet temperatures and thereby calibrate the nozzle 
thermocouple. 

In cases where the physical solubility of a gas in a liquid is limited, the slope of a plot of the 
inverse of the uptake coefficient (y) versus the square-root of the gas-droplet contact time yields the 
Henry’s law solubility constant and the intercept gives the gas-liquid mass accommodation 
coefficient. In Figure 4.2, yl is plotted against t '/2 for HC1 on 50 wt% H 2 SO 4 at four 
temperatures. Note that the slopes of the lines decrease as the temperature decreases, reflecting the 
increased solubility of a gas in a liquid at lower temperatures (slope «= H _1 , where H is the 
Henry’s law constant) . In addition, y(0) increases as the droplet temperature decreases, in 
agreement with our model for mass accommodation on aqueous surfaces [Watson et al., 1990] 
which predicts that the mass accommodation coefficient (a) increases as the surface temperature of 
the liquid decreases (we believe that the 251 K experiment is consistent with the other data within 
experimental error). In Figure 4.3, the logarithm of the effective Henry’s Law constant, H* (H* = 
H(l+Ka/[H+]), where K a is the HC1 dissociation constant) is plotted against T-f Note that the 
data point at 283 K represented by the open square was obtained from an earlier set of experiments 
[Watson et al., 1990]. In Figure 4.4, a plot of log(oc/(l-a)) vs. T-l is presented using data for 
HC1 uptake by 40 and 50 wt% sulfuric acid as well as data for HC1 uptake by water. 

Experiments on HC1 uptake by 40 wt 9c H 2 SO 4 droplets reveal that y is independent of HC1- 
droplet contact time (Figure 4.5). This is in contrast to the HC1 uptake data at 50 wt% (Figure 4.2) 
which vary markedly with contact time. This behavior is consistent with earlier experiments at 283 
K where it was determined that HC1 solubility in sulfuric acid is limited by the rate of the acid 
dissociation reaction HCl(aq) -> H+(aq) + Cl-(aq). The more concentrated the acid solution, the 



time {msec ) W2 

Figure 4.2. Plot of y 1 Versus t 1 / 2 for 50 wt% H 2 SO 4 
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Figure 4.5. Plot of y Versus t for 40 wt% H 2 SO 4 

less significant this reaction becomes. The flatness of the lines in Figure 4.5 indic ate s that the 
Henry’s law constants for 40 wt% H 2 SO 4 are too large to be determined from our data. In 
experiments on 70 wt% H 2 SO 4 droplets at 238 K, on the other hand, no HC1 uptake was 
detectable. Additional experiments are underway to complete the experimental matrix of 
temperature, wt% H 2 SO 4 , and contact time. These data will lead to the best values for the 
solubility of HC1 in H 2 SO 4 under stratospherically relevant conditions. 


STRATOSPHERIC AEROSOL THERMODYNAMICS 
Introduction 

Aerosols in the stratosphere are typically composed of aqueous H 2 SO 4 particles whose 
composition depends on the local temperature and water vapor pressures. This behavior is 
summarized in the phase diagram shown in Figure 4.6, which plots water vapor pressure versus 
temperature. The contour lines correspond to water vapor pressure curves of aqueous H 2 SO 4 
solutions (labeled by their H 2 SO 4 composition in wt %). TTie vertical scale on the right shows the 
corresponding altitude assuming a constant 5 ppm water vapor mixing ratio typical of stratospheric 
conditions, giving H 2 SO 4 solutions in the 60 to 80 wt % range at temperatures of 205 to 225K. 
These aerosol particles are typically < 0. 1 mm in diameter and believed to be in a supercooled 
liquid state. Total aerosol surface area peaks at about KF 9 cm 2 cm-3 in the 15-20 km altitude 
range, forming the Junge aerosol layer. This general picture has been quantitatively confirmed by 
extensive balloon measurements (Hofmann and Rosen, 1983, 1984), particularly after the El 
Chichon volcanic eruption which significantly enhanced aerosol densities (Hofmann, 1987). 
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Figure 4.6. Phase Diagram for H 2 O Vapor Pressure Versus T. Contour lines for Aqueous H 2 SO 4 
Solution (H 2 SO 4 Concentration in Weight %). 

Until recently it was believed that the only major exception to this behavior occurred near the 
poles where low temperatures in the winter led to the formation of PSCs containing ice particles 
(Stormer, 1929). One simply moved across the phase diagram to the circled x on the lower left of 
Figure 3.1. It was believed that as winter temperatures at high latitudes approached the frost point, 
PSCs grew by accretion of H 2 O on the background H 2 SO 4 nuclei (Steele et al., 1983). 

With the observation of large springtime ozone depletions over Antarctica, however, detailed 
analysis of observations from the SAM II instrument on Nimbus 7 (McCormick et al., 1982) 
provided evidence of significant aerosol enhancement at temperatures above the water frost point. 
Toon et al. (1986) suggested that HNO 3 might be a major component of these PSCs, and that 
significant quantities of HC1 might be incorporated as well. This led to attempts to construct phase 
diagrams for equilibria of H 2 O, HNO 3 and HC1 for stratospheric conditions (Wofsy et al., 1988; 
McEIroy et al., 1986). 

In the case of the binary H 2 O/HNO 3 system, these phase diagrams have been well reproduced 
by laboratory vapor pressure measurements originally reported by Hanson and Mauersberger 
(1988a, 1988b) and now confirmed in our laboratory (Worsnop et al., 1992). The experimentally 
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measured phase diagram is shown in Figure 4.7. The lines increasing diagonally to the left 
represent NAT vapor pressures obeying the relation 

p h 2 o 3 Phno 3 = k t 


The dashed lines represent previous measurements of Hanson and Mauersberger (1988a) while the 
points and solid lines are from our laboratory. The bolder diagonal lines increasing to the right 
show vapor pressures measured over samples of coexisting phases as labeled. The agreement 
between measurements from the two laboratories is excellent. 

The X in the center of Figure 4.7 corresponds to typical background HNO3 and H2O levels 
seen over Antarctica during field missions in August, 1987 (Fahey et al., 1989). That point 
intersects and NAT curve at 196 K, 6 K wanner than the 190 K frost point for ice found on the 
NAT/ice coexistence curve directly below it. Measurement of reduced stratospheric NO y densities 
in those 1987 field missions were indeed consistent with NAT condensation at higher temperatures 
relative to ice as predicted by the phase diagram in Figure 4.7. In fact, the boundary of the 
perturbed 



- 5 - 4 - 3 - 2 
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Figure 4.7. Phase Diagram for HNO3 and H 2 0 Partial Pressure as a Function of T for 
HN0 3 /H 2 0 Condensates. Points from Worsnop et al (1992); dashed lines from Hanson and 
Mauersberger (1988b). 
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a chemical region within the Antarctic vortex could by delineated by the NAT condensation 
conditions predicted by that phase diagram (Fahey et al., 1989). 

In contrast, detailed analysis of Arctic field observations in the stratosphere indicates there is a 
discrepancy between NAT phase diagram predictions and aerosol formation (Kawa et al., 1991a, 
1991b). Apparently aerosol free air is observed with HNO 3 vapor densities up to factor of 10 
higher than allowed by the NAT equilibrium relation (Kj). The authors hypothesize that this may 
reflect the degree of HNO 3 super saturation required to nucleate NAT on H 2 SO 4 aerosol particles. 
The apparent difference in temperature threshold for aerosol formation between the arctic and 
Antarctic can be explained by differences in temperature history of the sampled air particles, with 
those observed in the Antarctic spring having been colder for longer periods of time. A related 
question is whether the physical state of the background H 2 SO 4 aerosol plays a role. Under colder 
conditions, that aerosol may freeze, providing potentially better NAT nucleation sites relative to 
super-cooled liquid droplets (Dye et al., 1991). 

Our laboratory results support this supersaturation effect. The experimental method and these 
results are discussed in the following sections. 

Experimental Apparatus for Thermodynamic Measurements 

Our method for measuring vapor pressures in equilibrium with solid solutions of HNO 3 , HC1, 
H 2 O and H 2 SO 4 is based on infrared absorption using a tunable diode laser and a cryogenic 
temperature multiple-pass absorption cell. The apparatus shown in Figure 4.8 consists of a 4-inch 
diameter, high-vacuum, pyrex inner chamber which contains the multipass-cell mirrors. The solid 
sample is deposited at the bottom of the chamber on a cold spot which is ~6 mm in diameter. 

Cryogenic temperatures in the range 180 to 250 K are obtained by passing a regulated nitrogen 
gas flow through a heat exchanger immersed in a liquid nitrogen reservoir. This flow is then 
directed through a diffuser frit on the bottom of the surrounding stainless steel chamber. A heater 
on the diffuser frit regulates the temperature of a copper button which is epoxied onto the bottom of 
the pyrex inner chamber. A copper shield with a separate heater/sensor assembly is used to 
maintain the walls of the inner chamber and the mirrors at a temperature greater than the cold spot 
on the bottom. The entire assembly is enclosed in an evacuated stainless steel dewar for insulation. 

The tunable diode laser infrared absorption system uses a Laser Analytics cold head, current 
and temperature control modules. A laboratory microcomputer is used to scan the laser current and 
to signal average the detector output to obtain direct absorptions with a sensitivity of IO - 5 fractional 
absorption. The multipass mirrors in the cryogenic cell are of the off-axis resonator Herriott type 
and are set to obtain 40 passes with a base path length of 50 cm for a total pathlength of 2000 cm. 
The circular spot pattern is adjusted to minimize the amplitude of interference fringes arising from 
spill-over from adjacent spots (McManus and Kebabian, 1990). The gold coated pyrex mirrors are 
clamped concentrically to a 1 8 mm o.d. pyrex tube attached only to the top of the inner chamber. 
The clamping hardware is also gold coated so that only gold or pyrex surfaces are exposed to the 
gases during the vapor pressure measurements. Absorption line parameters and detection limits for 
selected gases are shown in Table 4-1. 
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Figure 4.8. Experimental Apparatus for Vapor Pressure Measurements 







Table 4-1 - Absorption Line Parameters and Detection Limits for Fugacity Measurements 


Line Position Line Strength^ Detection Limitb 

Molecule cm-l cm 2 molecule-1 cm-1 molecule cm-3 


HC1 

2821.570 

3.1 (-19) 

1.0(8) 

HBr 

2605.77 

1.0 (-19) 

2.0(8) 

HNO3 

1334.453 

3.4 (-20) 

3.1(8) 

H 2 0 2 

1302.506 

2.5 (-20) 

5.8(8) 

h 2 o 

1591.673 

1.7 (-21) 

1.5(10) 


a Line strengths are from AFGL line compilation (Rothman et al„ 1987) except for HNO3 
(May and Webster, 1987) and H2SO4 (Majkowski et al., 1978). 

b Based on a minimum detectable fractional absorption of AI/lo = 10-5 and a pathlength of 
2000 cm. 


All the exposed pyrex within the cells is covered with a hydrophobic Q>Fi3C2H4Si-(0)3- 
coating, which greatly reduces surface adsorption of H2O/HNO3 on the walls. This is critical to 
measurements of small trace gas concentrations, where it is important to eliminate wall effects 
which can dominate sample vapor pressures. A key aspect of the experimental design is the 
capability of independently controlling sample and absorption cell temperatures. Variation of the 
latter (controlled by the copper shield in Figure 4.8), while keeping sample temperature constant, 
can separate wall effects from sample vapor pressures. 

CURRENT RESULTS AND STRATOSPHERIC IMPLICATIONS 

HNO3 supersaturations are observed in our laboratory measurements as seen in Figure 4.7. 
For the 213 K data, the points lying directly above the vertical ice line (in lower right hand comer) 
correspond to HNO3 supersaturated over ice. We have observed HNO3 supersaturations > 5 that 
last 2-3 days over ice. Such observations are at least qualitatively consistent with the Arctic data 
discussed above. 

In fact, HNO3 supersaturation is repeatedly observed when an NAT free sample surface is 
prepared in our experiment. For example. Figures 4.9 and 4.10 show results of water addition to 
monohydrate (NAM) surfaces. In both figures the * points along an NAM line (above the 
NAM/NAT coexistence line) correspond to starting an experiment with a pure NAM surface. 
Adding H2O vapor to the sample cell then increases the H2O/HNO3 ratio of the sample surface, 
moving vapor pressures into the NAT portion of the phase diagram. In Figure 4.9, at 203 K, after 
H2O vapor addition to the cell, HNO3/H2O vapor pressures initially go to a metastable NAM/ice 
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Figure 4.9. HNO3/H2O Vapor Pressure at 203 K After Water Vapor Addition to NAM Sample 
(Labeled with *). Vapor pressures are close to a transient NAM/ice double poin t a fter 1 minute. 
Over the next hour (small circles) the vapor pressures relax to the NAT equilibrium line. Overnight 
(large circles) vapor pressures move toward NAM/NAT double point as the condensed surface 
sample equilibriates with the H2O/HNO3 ratio in the solid. 


double point, found by the extending NAM and ice vapor pressures (dashed lines). Over the next 
hour or so the vapor pressures "relax" to the NAT equilibrium curve as an NAT layer is nucleated. 

Figure 4. 10 shows the result of a similar experiment performed at 213 K. The data points 
show HNO3/H2O vapor pressures after successive small additions of H2O vapor. The points fall 
along a line of slope 1/2, which is consistent with nitric acid dihydrate (NAD) formation. For the 
sample shown in Figure 4. 10, further H2O addition led to "relaxation" to the equilibrium NAT line 
as in Figure 4.9. 

The behavior in Figures 4.9 and 4. 10 can be interpreted in terms of the degree of HNO3 
supersaturation required for NAT nucleation. In Figure 4.9, the NAM/ice double po int is 
supersaturated by over 100 in HNO3; thus, NAT nucleation began on a time scale of minutes. In 
Figure 4.10, condensation of a NAD surface layer occurred until HNO3 supersaturation of ~10 led 
to NAT nucleation. In both cases, metastable, transient phases (i.e., NAM/ice and NAD) were 
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Figure 4.10. HNO 3 /H 2 O Vapor Pressure at 213 K as a Function of Repeated Water Additions to 
NAM Sample (Labeled by *). Vapor pressures fall along a line of slope 1/2, corresponding to a 
dihydrate (NAD) phase. Further H 2 O vapor addition (not shown) led to NAT equilibrium (as in 
Figure 3.3), since HNO 3 supersaturation (~10) was sufficient to nucleate NAT. 


initially condensed on the surface. This behavior is typical of physical systems with significant 
free energy barriers to phase nucleation. 

These results can be summarized as follows: 

(1) At least under laboratory conditions, non-equilibrium behavior is typically observed for 
HNO 3 /H 2 O solid mixtures. For example, though not explicitly discussed above, the 
data indicates that surface and interior of solid samples are not equilibrated with the 
vapor conditions. Equilibration times of days for -5 mm thick samples correspond to 
solid diffusion coefficients on the order of KH 3 cm 2 sL 

( 2 ) Supersaturation factors of -5-10 in HNO 3 vapor are required for nucleation of NAT 
over ice, NAM and NAD surfaces. This is consistent with Arctic observations of 

apparent HNO 3 supersaturation. The key question is whether similar behavior is true 
for H 2 SO 4 surfaces (liquid or solid) that occur in the stratosphere. 




(3) Nitric acid dihydrate (NAD) is unstable with respect to NAT under stratospheric 
conditions. We have observed NAD at T = 214 K, but find no evidence that NAD is 
stable under stratospheric conditions. This indicates that the NAD crystalline phase 
observed via FTIR spectroscopy of thin films (Ritzhempt and Devlin, 1991; Koehler et 
al., 1991) probably does not occur in the stratosphere. 

(4) Under stratospheric conditions the equilibrium phase of the HNO3/H2O system is pure 
NAT (or NAT/ice mixtures for T < 190 K), consistent with the thermodynamics 
implied by the Kj vapor pressure relation. Based on our observation of slow solid 
phase diffusion, we believe earlier reports of "water rich" NAT samples (Hanson and 
Mauersberger, 1988b) came from non-equilibrated solid samples. 


RELEVANCE OF EXPERIMENTAL RESULTS TO 
STRATOSPHERIC FLIGHT 

The overall goal for the HSRP program is to integrate variations in aerosol composition and 
distribution with the heterogeneous reaction chemistry represented by R1-R6. 

The growing diversity of aerosol observations and apparent complexity of NAT nucleation is 
very relevant to the modeling of potential chemical perturbations due to aircraft emissions. Those 
emissions are projected to occur in the center of the background sulfate layer. Especially at lower 
temperatures at lower altitudes (< 20 km) near the tropopause and at higher latitudes in die winter, 
emission of H2O vapor and condensation nuclei could significantly enhance aerosol formation, 
since background conditions there are already very close to NAT condensation predictions (Hamill 
and Fiocco, 1988). Combined with reports of possible observation of local ozone perturbations 
due to direct interaction with Arctic and Antartic aerosols (Hofmann, 1989; Hofmann et al., 1989; 
Fiocco et al., 1989), it is important to quantitatively understand the vapor pressure diagrams for 
prediction of aerosol formation and growth. 

Finally, it is also important to extend these vapor pressure measurements to trace vapor 
pressures over aerosol solutions containing HCI/H2O/HNO3/H2SO4. The measurements of 
Hanson and Mauersberger (1988c, 1989) of HC1 solubility in ice and NAT solutions are indicative 
of the difficult nature of these experiments. These measurements are directly relevant to 
heterogeneous reactions such as R5-R6. To some extent, at higher temperatures these solubility 
measurements can provide important overlap with uptake experiments described in the previous 
section. 
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of the condensed solid water (ice) particles and on exhaust soot particle surfaces. 
Dispersion and reaction of chemical constituents in the far wake are studied with a 
Lagrangian air parcel model, in conjunction with a radiation code to calculate the net 
heating/cooling. Laboratory measurements of heterogeneous chemistry of aqueous 
sulfuric acid and nitric acid hydrates are also described. Results include the 
solubility of HC1 in sulfuric acid which is a key parameter for modeling stratospheric 
processes. We also report initial results for condensation of nitric acid trihydrate 
from gas phase H 2 O and HNO^. 
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